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FRACTIONAL CAHN-HILLIARD, ALLEN-CAHN AND POROUS MEDIUM 

EQUATIONS 

GORO AKAGI, GIULIO SCHIMPERNA, AND ANTONIO SEGATTI 


Abstract. We introduce a fractional variant of the Cahn-Hilliard equation settled in a bounded 
domain D C M. N and complemented with homogeneous Dirichlet boundary conditions of solid type 
(i.e., imposed in the whole of IP ;V \ O). After setting a proper functional framework, we prove 
existence and uniqueness of weak solutions to the related initial-boundary value problem. Then, we 
investigate some significant singular limits obtained as the order of either of the fractional Laplacians 
appearing in the equation is let tend to 0. In particular, we can rigorously prove that the fractional 
Allen-Cahn, fractional porous medium, and fractional fast-diffusion equations can be obtained in 
the limit. Finally, in the last part of the paper, we discuss existence and qualitative properties of 
stationary solutions of our problem and of its singular limits. 


1. Introduction 

Let fl be a smooth bounded domain in M. N . For s,cr £ (0,1), we consider the following class 
of initial and boundary value problems: 


(1) 

dtu + (—A ) s w = 0 

in x (0, -|-oo), 

(2) 

w = (—A ) a u + W'(u) 

in x (0, + 00 ), 

(3) 

u(x, 0) = uo(x) 

in fl, 

(4) 

u = w = 0 

in \ n. 


The above system constitutes a natural generalization of the well-known and extensively studied Cahn- 
Hilliard equation , Allen-Cahn equation , and porous medium equation. More precisely, when s = a = 1, 
system (l)-(2) reduces to the Cahn-Hilliard equation (cf. [12]), when s = 0 and cr = 1 we have the 
Allen-Cahn equation (cf. [3]), and when s = a = 0, (l)-(2) turns to an ODE. The relations between 
the above system and the porous medium equation will be outlined and made rigorous later on. 

The function W in (2) represents a configuration potential which may have two (or more) 
wells. The general structure of W is given by 

(5) W{v) = j3(v) — ^u 2 , so that W'(y) = /3(v) — Xv, where /3 = /?'. 

Here j3 is a smooth and convex function and A > 0 is a constant. Hence, if A > 0, then W may be 
nonconvex. In the phase-transition literature, the wells of W correspond to energy minima (attained 
at pure phases or configurations). In view of the variational structure of our system, it is convenient 
to keep the same interpretation also in the present case. Actually, several types of significant choices 
have been proposed for W, also including cases where j3 is nonsmooth or even singular (like the so- 
called logarithmic potential $(v) = (1 — v) log(l — u) + (1 + v) log(l + v), v £ (—1,1)). For the sake 
of simplicity, in this paper we will just consider the case given by 

(6) f3(v) = -\v\ p , p£( l,oo), so that f3(v) = |f| p-1 sign v. 

Moreover, we will set A = 1. This choice, up to an additive constant, includes the standard double-well 
potential, namely W(v ) = j{v 2 — l) 2 , widely used in the literature. Mathematically speaking, the 
cases p > 2 and p < 2 enjoy rather different features. Indeed, while in the former case W ( v ) is coercive 
as |u| oo, in the latter situation W is unbounded from below. Hence, we can expect (and, indeed, 
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we will prove) different behaviors of the solutions, especially for large values of the time variable. Note 
that the case p = 2 corresponds in fact to the linear problem and is not considered here. 

A further important feature of our problem is the occurrence of solid boundary conditions 
of homogeneous Dirichlet type, stated by (4). Namely, the values of u and w are prescribed in the 
whole complement of f2, not only on the boundary. Of course, this assumption is strongly related to 
the nonlocal character of fractional Laplacians. Indeed, it is worth observing from the very beginning 
that, though the values of (— A) s w and of (—A ) a u at any point x £ fl depend also on the values of 
u and w outside Q (which are set to be 0 by the boundary conditions), we prescribe the validity of 
(1) and (2) only at the points x £ fl. For x £ 1^ \ H, (l)-(2) need not to be satisfied (and, indeed, 
there is no reason why they should). This observation will be further clarified in Section 2, where the 
appropriate concept of weak solution is introduced. Correspondingly, we will also recall the precise 
definition of (—A) r , r £ (0,1), both in the strong and in the weak (variational) form, the latter being 
the more appropriate one for the analysis of our problem. 

The study of system (l)-(4) is motivated both from the point of applications and due to its 
mere mathematical interest. Under the first perspective, it is worth recalling that the Allen-Cahn and 
Cahn-Hilliard equations, in their standard formulation (i.e., with the usual Laplace operators), play 
a central role in materials science. Indeed, they commonly occur in mathematical models for phase- 
transition or separation, viscoelasticity, damaging, complex fluids, and whenever diffuse interfaces 
appear (see, e.g., [13, 28] for a comprehensive bibliography). A reason for considering a fractional 
version of the Cahn-Hilliard equation can be provided by observing that, in the original formulation 
of the physical model [12], the Laplace operator in (2) was actually replaced by a spatial convolution 
term, aimed at describing long-range interactions among particles. It was only in the subsequent 
mathematical literature that, mainly for analytical reasons, this nonlocal term has been substituted 
with the term —A u. Under this perspective, the use of the fractional Laplacian (—A) CT (which, at 
least for smooth functions, may be represented exactly by a convolution integral) appears to be more 
adherent to the physical setting. It is worth remarking that the study of fractional (or, more generally, 
nonlocal) PDE’s is a lively research topic, both from the point of view of mathematical theory and 
in relation with the many real-world applications. Among these, we mention obstacle problems [10], 
finance [16], quasi-geostrophic flows [11, 15], anomalous diffusion [27, 37, 39]. A more comprehensive 
list of references is provided in the survey [18]. 

In the recent literature, a relevant number of works have been devoted to the analysis of 
nonlocal Allen-Cahn and Cahn-Hilliard models. Here we quote, among others, [4, 20, 29] (see also 
[14] for an application to complex fluids). Actually, in [4, 14, 20], the term (—A ) a u in (2) is replaced 
by 

(7) J[u\ = a(-)u - j * u, a(x) = / j(x-y)dy, (j*u)(x)= / j(x - y)u(y) dy, 

Jn Jn 

where the convolution kernel j enjoys suitable, and rather strong, regularity properties. For instance, 
in [4], j is assumed to lie in C 2+ “(K Ar ) for some a > 0, while in [14] j is taken in IU 1,1 (]R Ar ), and in 
[20] a similar condition is required. In all cases, both integrals in the right hand side of (7) are finite 
for (almost) every x £ S2, whenever u lies, say, in L 1 (U) (more comments on this point will be given 
in the Appendix). 

On the other hand, the kernel K r generating the fractional Laplacian (—A) r (cf. (10) below) 
is not even summable. Hence, the corresponding convolution integral needs to be intended in the 
principal value sense even for smooth functions u (cf. (11)). In addition to that, if u is not regular 
(for instance if it just lies in some L p -space), then the “pointwise” expression (10) makes no sense at 
all and a variational definition of (—A ) r u is required. In this sense, the fractional Laplacian gives rise 
to a much stronger singularity with respect to those considered in [4, 14, 20]. Up to our knowledge, 
the only papers dealing with a true fractional Allen-Cahn model are [29] and [19], where, however, 
the analysis is restricted to the spatial one-dimensional case and mostly concentrated on other aspects 
rather than weak solvability and regularity of solutions. 

More recently, Abels, Bosia and Grasselli in [2] analyzed a variant of problem (l)-(2), where 
the diffusion operator in (the analogue of) (1) is the standard Laplacian, while the diffusion operator 
in (the analogue of) (2) is the so-called regional fractional Laplacian. For its definition we refer the 
reader to the Appendix; here it is just worth mentioning that its properties are slightly different 
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compared to those of the operator considered in this paper, the main point regarding the boundary 
conditions. Actually, the operator of [2] can be seen as a fractional power of the Neumann Laplacian 
on Q (in particular, this implies conservation of mass, which does not hold here in view of the Dirichlct 
condition (4)). The authors of [2] prove existence and uniqueness of weak solutions in the case when 
the function W may have a singular character (as happens for the logarithmic potential mentioned 
before). Moreover they characterize the long-time behavior of solution trajectories proving existence 
of the global attractor for the dynamical process associated to the system. 

The first aim of this paper is to prove existence and uniqueness of solutions to a weak formula¬ 
tion of problem (l)-(4). Hence, differently from [2], we will consider a fractional dynamics both in (1) 
and in (2). This program requires, at first, to set a proper functional framework. Our approach basi¬ 
cally follows (and complements) the perspective given in [33], where a weak version of the fractional 
operators ruling the system is defined. Actually, a variational expression of the Dirichlet fractional 
Laplacian can be given at least in two (equivalent) ways: one can either work with functions defined 
on Q and implicitly extend them to 0 outside fl when computing fractional Laplacians (which depend 
on values taken on the whole of R^), or one may use spaces of functions defined on the whole space 
but constrained to be identically equal to zero outside H. Generally, we shall work within the latter 
framework. This choice permits us to address the problem, for all fixed s,er £ (0,1), in the usual 
Hilbert setting, very similarly to what happens for the standard Cahn-Hilliard model. In particular, 
we can prove existence by means of a classical time-discretization scheme. Compactness and dual¬ 
ity arguments are then exploited in order to pass to the limit in the discretization; in this way we 
can avoid any reference to finer regularity properties of solutions to fractional elliptic and parabolic 
problems, which may involve rather delicate issues. Uniqueness also follows from a simple contraction 
principle. 

As anticipated, after establishing well-posedness of the model, we will turn our attention to 
further properties of solutions. As a first issue, we will let the “order” er of the fractional Laplacian 
in (2) go to 0. This singular limit is motivated by noting that, as er \ 0, (— A) a u —> u in a suitable 
sense. As a consequence we obtain, at least formally, 

w = (-A Yu + W'(u) = (-A Yu + /?(u) - u p (u). 

Hence, on account of (1) and recalling assumption (6), one expects to get in the limit the equation 

(8) d t u + (—A) s (|it| p_1 signu) = 0 in fl x (0,+oo), 

(9) u = 0 in K.^ \ fl. 

This corresponds, for p > 2, to the so-called fractional porous medium equation , recently addressed 
and studied in a number of contributions (see, among others, [17], [5], [6] and [7]). 

Actually, taking a family { (u a , w a )} of solutions to our problem, and letting er \ 0, we can 
rigorously prove that, up to extraction of a subsequence, u a tends to a limit function u satisfying (8). 
Our result holds, under natural assumptions on the initial data, both for p > 2 and for p £ (1,2). 
The proof is, however, not straightforward and relies on some fine properties of first eigenvalues of 
fractional elliptic Dirichlet problems, which, to the best of our knowledge, are new (see Proposition 
2.2). The case p £ (1,2) is a bit more involved due to the lack of coercivity of the energy functional. 
Indeed, in that case we need to modify a bit the energy (see Theorem 3 below) in such a way to get an 
estimate for u a uniform as a 0. However, this modification does not affect the limit equation (8), 
which corresponds in this case to the fractional fast-diffusion equation studied, e.g., in [24]. 

On the other hand, it is also natural to investigate what happens as one lets s 0. In that 
case, one expects that 

(-A ) s w s -^4 w. 

In other words, the limit gives rise to the fractional Allen-Cahn equation 

d t u + (—A ) a u + W'(u) = 0 inflx (0, +oo), 
it = 0 mR N \Q, 

both for p > 2 and for 1 < p < 2. It is worth noting that the limit s \ 0 is considerably simpler than 
the limit a \ 0 since keeping er fixed maintains some additional space compactness, which reveals 
to be helpful for the purposes of obtaining a strong convergence for u and identifying the nonlinear 
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terms in the limit. As a result of these two procedures, we can see problem (l)-(4) as a bridge between 
the usual Cahn-Hilliard equation (given by (<r, s) = (1,1)) and the fractional (or also non-fractional) 
porous medium and Allen-Cahn equations. 

The last part of the paper is devoted to the proof of some results related to stationary solutions 
to problem (l)-(4) (for fixed s,cr £ (0,1)). As expected in view of the nonconvex character of W, we 
can show that nontrivial stationary states exist if and only if the first eigenvalue Ai(cx) of (—A) CT is 
strictly smaller than 1. Indeed, whenever Ai(cr) > 1, the coercivity given by (—A)' 7 compensates the 
nonconvexity of W'(u) = fi(u) — u, exactly as happens for the standard Laplacian. Of course, the 
properties of the stationary states play an important role for what concerns the long-time behavior of 
solutions to the evolutionary system. We plan to address this issue, and, particularly, to investigate 
the properties of ui-limit sets , in a forthcoming paper. 

The plan of the paper is as follows. A survey on some basic definitions and tools related 
to the fractional Laplacian, as well as some useful lemmas, are presented in the next Section 2. In 
Section 3, we introduce our assumptions and state our main results. The proofs of existence and 
regularity properties of solutions are carried out in Section 4, while the convergence to the fractional 
porous medium and Allen-Cahn equations is analyzed in Section 5. Our last results regarding the 
properties of stationary states are presented in Section 6. Finally, in the Appendix we provide some 
more comments on the relations occurring between the problem analyzed here and other nonlocal 
models of Allen-Cahn or Cahn-Hilliard type studied in the literature. 


2. Notations and Preliminaries 


2.1. The fractional Laplacian. We introduce here the standard (strong) form of the fractional 
Laplace operator (—A) r in the whole space R N . The reader may refer, e.g., to [18] for additional 
details. Given r £ (0,1), for u in the Schwartz class S(R N ) of the rapidly decaying functions at 
infinity, (—A ) r u is defined as 

u(x) - u(y) 


( 10 ) 


(—A ) r u(x) := C(r,N)p.v. [ 

Jr 1 


■dy, 


i rm \ x - y\ N+2r 

where the notation p.v. means that the integral is taken in the Cauchy principal value sense, namely 


( 11 ) 


p.v. 


u(x) - u(y) 
-nv+27 d v = 


u{x) - u(y) 


d y. 


jk- \x - y\ N+2r Ko J R N\ B (x,e) \x - y\ N+2r 
The exact value and the asymptotics with respect to r of the normalizing constant C(r, N ) are crucial 
for our purposes. To this end, we recall that C(r,N) = (f RN d^) -1 and that (see, e.g., [18, 

Corollary 4.2]) 

, 19 x r C(r,N) 2 

(12) !\or(l-r) 

On the other hand, since the dependence of C(r, N ) with respect to the space dimension N is not the 
major issue for this paper, we will always write C(r) for C(r,N) in what follows. For any r £ (0,1) 
and for any x,y £ R w we will also use the shorthand notation 

K r (x — y) = \x-y\~ N ~ 2r 

to denote the singular kernel appearing in the definition of (—A) r . A second, albeit equivalent, 
definition can be given using the Fourier transform. Indeed, (—A) r can be introduced as the pseudo- 
differential operator of symbol |£| 2r , namely 

(13) (-Arn = ^- 1 (|e| 2r j?(n)), V«g 5(R jv ). 

We denote by $(v) (or by v) the Fourier transform of v. 


2.2. Fractional Sobolev spaces. In this subsection, we shall deal with fractional Sobolev spaces. 
We refer the reader to, e.g., [25] and [1] for further details. 

For r £ R, the fractional Sobolev space H r ( l w ) is defined by 

H r (R N ) := {n £ S(R N )': (1 + |e| 2 ) r/2 i)(C) £ L 2 (Rf )} , 








FRACTIONAL CAHN-HILLIARD, ALLEN-CAHN AND POROUS MEDIUM EQUATIONS 


5 


where S(R N )' stands for the dual space of the Schwartz class S(R N ) and L 2 { R^) is the space of 
square-integrable functions with respect to the variables £ £ R^, equipped with the norm 


IMh= (i + |£| 2 ) r/2 ®(0 


L 2 ( Rf) 


for v £ H r {R n ). 


In particular, for r £ (0,1), we can equivalently write 

H r {R n ) := {v £ L 2 (R n ) : (x, y ) K r (\x - y\)\v(x) - v(y )| 2 £ L\ R 2W )} , 
endowed with the (equivalent) norm 


(14) 


MIlu-(rn) IM| 2 2(R«) ■+-f° r v € H r (R N ). 


Here, [•]#•»• denotes the so-called Gagliardo-seminorm 


lH r (R N ) 


// K r (x - y)\v{x) - v(y)\ 2 dxdy. 

J J R2JV 


Furthermore, for r £ R, H r (R N ) is an intermediate space between H m ( R") and L 2 (R 7V ), that is, 


H r 


r,N 


) = [FT m (R JV ),L 2 (R 7V )] 


for any 9 £ (0,1) and m £ Z satisfying r = (1 — 8)m. 

Let U be a bounded domain of R N with smooth boundary dfl. For r £ R, the fractional 
Sobolev space H r (Q) may be analogously defined as 

H r m := [H m (n),L 2 (n)] g 

for any 6 £ (0,1) and m £ N satisfying (1 — 9)m = r. Moreover, for r £ (0,1), one may use an 
alternative definition, 

H r (fl) := {u £ L 2 (VL ): (x, y) i —> K r (\x — y|)|zi(a;) — u(y )| 2 £ L 1 (H x H)} 
with the intrinsic norm 

C(r) / 


\\ v \\h-(Q) : =\\ v \\l^n) + ^ip- I[ Kr{\x - y\)\v(x) - v(y)\ 2 dxdy for v€H r {n). 


Sobolev embeddings and inequalities also hold for fractional Sobolev spaces. Thus, H r ( fl) is con¬ 
tinuously (resp., compactly) embedded in L q (VL ), provided that 1 < q < 2* := 2N/(N — 2r) (resp., 
1 < q < 2*) and 2r < N. Finally, for each r > 0, Hq(Q) is defined as the closure of in U r (f2). 

In case r < 1/2, the space Hq(£1) coincides with H r ( fl); in case r > 1/2, is strictly contained 

in H r (Q) (see, e.g., [25, p. 55, Theorem 11.1]). 


2.3. The functional framework. It is apparent from (13) that, for v £ ^(R^), (— A) r v does not 
necessarily belong to S(R N ) (being r < 1, the symbol |£| 2r introduces a singularity in the origin in 
its Fourier transform). Moreover, even for v with compact support, (—A)T generally does not have 
compact support due to the non locality of the operator. In addition to this, the above definition could 
make no sense when non-smooth functions are involved. Thus, it will be important for us to extend 
the definition of the fractional Laplacian to a more general setting. This will be accomplished by using 
the theory of distributions together with some tools of convex analysis. The framework we are going 
to fix will permit us to use variational and energy techniques in order to address our problem. As 
already observed in the introduction, although equations (l)-(2) are settled only in H, the behavior 
of (—A ) a u and (—A ) s w depends on the interplay between the values of u and of w inside and outside 
Q. Proceeding along the lines of [33], we can then introduce some functional spaces. 

Firstly, we set 

(15) H 0 := {i) £ L 2 (R n ) : v = 0 a.e. in R N \ H}. 

The space L 2 (R N ) (hence its closed subspace Ho) is endowed with its standard scalar product, 

(u,v) := / u(x)v(x)dx for u, v £ L 2 (R N ). 

Jr n 
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Of course, in the closed subspace H 0 taking the scalar product of L 2 ( Cl) and the associated norm would 
make no difference. Hence, L 2 (Cl) can be identihed with Hq by zero extension outside Cl. Furthermore, 
we set 

L p 0 (R n ):={vGL p (R n ):v = 0 a.e. in R w \Cl} with || - IU S (a^) := || - 

^i(R*) : = {0 g C°°(R n ) : <t >|n £ C^°(Cl) and 0 = 0 in R N \ (1 } , 

which can be identified with L p {Cl) and Cq°(Q), respectively. Then Lq(R n ) is reflexive for p £ (1, oo), 
since L P (R N ) is closed in L P (R N ); moreover, Lq(R n ) is separable for p £ [l,oo). Furthermore, 
@n(R N ) is dense in L^R^), provided that 1 < p < oo. 

As in [33], we set Q := R 2Ar \ ((R w \ fl) x (R N \ Cl)) and denote by X r _o, r £ (0,1), the space 

(16) X rfi := |u £ H 0 : (x,y) t—>■ (v(x) - v(y))y/K r (x - y) £ T 2 (Q)j. 

Actually, X r $ can be endowed with the scalar product 

C(r) [ f 

(17) (v,z) Xr0 :=(v,z) + —^-jj K r (x — y)(v(x) — v(y))(z(x) — z{y)) dxdy for v,zGX r>0 
and the associated norm 

(18) IMI 2 := IMlia(n) + -ip-JJ Q K r(x-y)\v{x)-v(y)\ 2 dxdy for v £ X r>Q , 

where C{r ) is as in (10). Then, it is easy to check that X t q is a Hilbert space (i.e., the above norm 
is complete). Note that X r p could be also presented in a more familiar form, 

(19) X r ,o = {u £ H r (R N ) such that v = 0 a.e. in R w \ fl}. 

There holds the Poincare-type inequality, 

(20) M 2 L 2 {n) <c P (r)^-[v} 2 Hr{RN) for all v £ X r> o 

for some constant cp(r) depending only on r, N and the diameter of Cl. Indeed, one can take R > 0 
such that Cl is included in the open ball Bp of radius R centered at the origin. Then, using the 
definition of the Gagliardo-seminorm, we see that 

Mjp-(rJV) > / / K r {\x-y\)v{x) 2 dxdy 

J n c Ja 

> b K r (\x - y\) dyj v(x) 2 dx > (2^2)^+^ for a11 V £ * r ’ 0 ’ 

where stands for the complement of Cl and \Bp + \ \ fl| denotes the Lebesgue measure of the set 
Br + i\CI. Note that > \Bp + i\Bp\ >0. Thus the quotient M#r( R iv)/|M| 2 2 (Q) is bounded 

from below for v £ X t q, whence (20) follows. Hence, by (20), the norm on X t q given by 

(21) \\ v \\ 2 x r , 0 '■= JJ^Kr(x-y)\v(x)-v(y)\ 2 dxdy for v £ X r>0 

is equivalent to that defined in (18). From now on, we will fix (21) to be the norm in X r ^. 

Now let us introduce the dual spaces H' Q and X' rQ of Hq and X r p, respectively. In particular, 
we are here concerned with the meaning of the equality in the dual space , 

f = 9 in H' q for /,<? £ H' 0 , 

which actually means 

(22) (f,v ) H o = (g,v) H o for all v £ H 0 . 

By the Riesz representation theorem, one can uniquely take Uf,u g £ Hq such that {f,v)n 0 = ( u/,v ) 
and ( g,v)n 0 = {u g ,v) for all v £ H 0 . Hence (22) yields Uf =u g in Cl (hence, over R N ). On the other 
hand, one may generate h £ Hq from a function Wh £ L 2 (R N ) which might not vanish outside Cl by 
setting 
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Then h coincides with / in Hq, provided that u>h = Uf in Cl. In other words, even if / and g have 
pointwise representations Wf,w g £ L 2 (R JV ), respectively, as in (23), the relation / = g in H' 0 ensures 
that Wf(x) = w g (x) for a.e. x £ fi only, and it does not guarantee the coincidence of w/ and w g 
outside B. This observation is also extended to the relation f = g in X^ Q . 

From now on, we identify the Hilbert space Hq with its dual space Hq by means of the scalar 
product (•, •). More precisely, we shall identify / £ Hq with its unique representation it/ £ Hq by the 
Riesz representation theorem. Hence we shall write / = g in Hq for claiming the equality of /, g £ Hq 
(i.e., (22)) as well. Here we should emphasize again that whenever are given representatives Wf and 
w g of /, g £ Hq , respectively, the identification implies that wf = w g only in Q. Then since X r ^ can 
be seen as a dense subspace of Hq, one may consider the Hilbert triple, 

(24) Xr t0 ^Hq-H'q^t X' q, 

with compact and densely defined canonical injections. This relation will play a crucial role throughout 
this paper. 

On the other hand, for r £ (0,1), the extension operator of u £ X r $ to 0 outside Q is a 
continuous mapping of H r {Vl) —> H r (R N ). In particular (see [25, Theorem 11.4, Chapter 1]), if 
r £ (1/2,1), the functions in X r ^ are equal to zero, in the sense of traces, on dfl. Hence, X t q can be 

identified with Hq(Q) in that case, whereas X r Q ~ Hq(JI) = H r (Q) for r £ (0,1/2). Finally, in the 

1 /2 

limit case r = 1/2, it turns out that X\/ 2 ,o — H 0 q (Q) (again, see [25] for more details). 

Based on this functional framework, we can introduce, for r £ (0,1), the weak form 2l r of the 
fractional Laplacian (—A) 7 '. More precisely, the operator 2l r : X r $ —> X{, 0 is defined by 

(25) (21 r v,(j>) :=^p- [ [ K r (x-y)(v(x) - v(y))(<j)(x) - c/>(y))dxdy, for all v, <j> £ X r<0 , 

where the integral over M. 2N can be equivalently replaced with an integral over Q. Note that, as soon 
as v, 4> £ X r Q, the integral in the right hand side is finite. Note also that (25) can be understood as an 
integration by parts formula, at least when v, (j) are sufficiently regular. Indeed, to see this, we define 


(26) 


(-A)Mr) := C(r) [ K r>e (x - y){u{x) - u(y)) d y, 
Jr* 


where K r e := K r { 1 — Xb(o,e )) an< i Xb(o,e) denotes the characteristic function of the ball B{ 0,e) in 
R^. Then, by symmetry of K r e , 


(27) 


{-A) r s v(x)<l)(x)dx = C(r) [[ K r , s (x - y)(v(x) - v(y))<j)(x) dx dy 

J Jr 2N 


Jr n 

— C(r) [[ K r £ (z)(v(x) — v(x + z))(f>(x) dx dz 

J Jr 2N 

= C(r) [j K r , e (z)(v(x) — v(x — z))(j>(x) dx dz 

J J R2JV 


^ f f K r e (z)(v(x) — v(x — z)){(j>{x) — 4>{x — z)) dx dz 

J J R 2JV 

/ / Kr tS (x - y){v(x) - v{y)){(j)(x) - <j>(y))dxdy. 

J Jr 2N 


2 

C(r) 


Now, letting er \ 0, the right hand side converges to (2l r e, <j>), provided that v, <fi lie in X r 0 (and, hence, 
a fortiori if v, <f> are smooth functions). On the other hand, whenever v is so smooth that (—A) r v (i.e., 
the “strong” fractional Laplacian of v defined in (10)) is represented by, say, an L 2 -function, then 
the left hand side of (27) converges to J rN (— A) r v(x)cj)(x) dx. Hence, 2t r can indeed be seen as an 
extension of (—A) r to less regular function. Moreover, formula (25) can be also expressed in terms of 
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Fourier transform. Actually, for v, </> £ X r q, thanks also to Fubini’s theorem, we have 

(2l r u, cj>) = f ( f K r (z)(v(y + z)~ v(y))((/>(y + z) - cj>(y)) d?/) dz 

z J i« \Jrn J 

= ^T f HOW) ( [ K r {z) |e** - 1| 2 d z) d£ 
z J i« \Jr n J 

( 28 ) = / Kr«(£)jiFI(£)d& 

since J Rjv K r (z)\e l ^' z — 1\ 2 dz = |£| 2r (see [18]). In particular, using (13), we get 

(29) <2U,^}=/ \Z\ r mWfa)<%= f (-A) r / 2 v(x){-A) r/2 <t>(x)dx = (X r ct>,v), 

Jr n Jr n 

which could serve as an alternative, albeit equivalent, definition of the weak fractional Laplacian 2 l r . 

A weak form of the fractional Laplacian can be introduced also for the whole space case, 
fl = l w . Indeed, for r £ (0,1), we can set (cf. (16)) 

(30) := {v £ L 2 ^) : (x,y) ^ (v(i) - v(y))y/K r (x-y) £ L 2 (R 2Ar )} = H r (R N ). 

Then, the previous discussion extends, with minor modifications, to the space X r . Moreover, one 
can correspondingly consider the Hilbert triplet (X r , H, Xf). With a small abuse of notation we will 
indicate with the same symbol 2l r the weak form of the fractional Laplacian as an operator from X r 
to X' r defined by (25) with X r $ replaced by X r . 

The relations between the Gagliardo-seminorm and the Fourier-transform definition of the 
fractional Laplacian are also clarified by the following property (cf. [18, Propositions 3.4 & 3.6]): 

(31) = 1124 / 2 ^ 11 !^) = for v £ H r (R N ) and r £ (0,1). 

In particular, the norm in H r (JSL N ) can be equivalently expressed as 

(32) IMIif(]R JV ) = II^IIl 2 (R jv ) + l|2tr/2^||i2(]RJV) = IMIl 2 (R n ) + ||l£l 

It is worth noting that there is another possible approach for dealing with fractional Laplacians on 
bounded domains. Indeed, one may define the operator, called spectral fractional Laplacian , as 

+oo 

(33) (—A n ) r f(x) := ^ A T j fj4> j {x), x £ fl, 

.7=1 

where Ay > 0,j = 1,2,..., are the eigenvalues of the Dirichlet Laplacian on f2, <f>j are the corresponding 
normalized eigenfunctions, and 

fj ■= / f(x)<t>j(x)dx, with \\(f)j\\ L 2 ^) = 1- 

Jo. 

As observed in [34], the spectral fractional Laplacian (—An) r is a different operator with respect 
to the operator considered in this paper. We refer to [9] and [7] and to the references therein for 
the functional framework related to (33) and for the analysis of some differential problems involving 

(-Aar- 


2.4. An L 2 -framework for fractional Laplacians. The weak fractional Laplacian 21, can be also 
interpreted in the framework of convex analysis. Actually, for r £ (0,1), we can introduce the 
functional G r : Hq —> [0, +oo] given by 


(34) 


G r (v) 


JJ Q Kr ( x - y )H x ) 

+00 


v{y)\ 2 dx dy 


if v £ X rfi , 
otherwise. 


Then, it is obvious that G r is a lower semicontinuous (in Hq), convex functional. Moreover, given 
u £ D(G r ) := X t q (the effective domain of G r ), it is clear that, for all v £ X r ^, 

.. G r (u + tv) — G r (u) 

hm--- 

t —>-0 


(35) 


t 


(21 r u,v) Xr 0 . 
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On the other hand, if £ £ H 0 belongs to dG r (u ), where the subdifferential dG r is defined by 

dG r [w) := {£ £ H 0 : G r (v) — G r (w) > (£, u — w) for all v £ H 0 } for w £ D(G r ) 

with domain D(dG r ) := {tu £ D(G r ): dG r (w) ^ 0}, then one can easily check that, again for all 
n £ X r< o> 


(36) 


G r {u + tv) - G r {u) , , 

hm-= (£, v). 

t—>0 t 


Combining (35) and (36), for u £ D(dG r ) and £ £ dG r (u), we obtain 


(37) (£, v) = {Ql r u, v) Xr o for all v £ 7f ri0 , 

whence dG r (u ) = {£} and £ is the realization in Hq of 2t r u. Then is unbounded linear in Hq. 
From now on, dG r will be denoted, with a small abuse of notation, by 2t r . Then D( 2l r ) means D(dG r ); 
moreover any u £ D( 2l r ) can be seen as the solution to the elliptic problem 


(38) 


2 l r '« = / in i7 0 


for some / £ 77o- Here we note that (38) does not mean that 21 r u vanishes outside Q (see §2.3). Due to 
the lack of regularity results for fractional elliptic problems such as (38), it is a nontrivial issue to give 
a precise characterization to the domain D( 2l r ). Up to our knowledge, the sharpest results available 
in literature are due to Ros-Oton and Serra, who proved in [31, Prop. 1.4 (ii)-(iii)] the following: 

Proposition 2.1. Let LI C be a bounded C 1,1 -domain. Then, if r £ (0, -j) D (0,1), the solution 
u to (38) satisfies 

(39) IMU?(fi) < C\\f\\ Ho for q = ^ 
while for r £ (-j, 1) fl (0,1) we have 

(40) IMIc<* ( n) < C'll/lko f° r a = min|r,2r - y| . 

In both cases the constant C > 0 depends only on r, |H|, and q (or a). 


The operator 2l r also enjoys a useful monotonicity property. Indeed, let /3 : R —> R be any 
smooth monotone function such that /3(0) = 0 and let $ : R —> [0, +oo] be the (convex) function such 
that /3(0) =0 and $' = ff Then, setting for instance 

! P(r) if |r| < e -1 , 

/9(e _1 ) if r > e _1 , 

/3(—e _1 ) if r < —e _1 , 

then, for fixed e > 0, /3 e is bounded, Lipschitz continuous, and monotone. Hence, it is immediate to 
check that, if v £ X r ,a, then f3 e (v) £ <T r ,o for each e > 0. Hence, by monotonicity, 

(42) (K r v,p e (v))x r , 0 = [I K r (x - y)(v(x) - v(y))(ff(v(x)) - P E (v(y))) dxdy > 0. 

z J ii 2N 

Moreover, it is clear that, for any n,r 2 £ K, the product (/3 e (r\) — /3 e (r2))(ri — 7 - 2 ) increases as £ 
decreases. Hence, letting e \ 0, we can apply the monotone convergence theorem to (42). This gives 

(43) —tP f [ K r (x-y)(v(x) - v(y))(ff(v(x)) -/3 e (v(y))) dx dy 

1 J J R2JV 

-t —P [[ K r (x - y)(v(x) - v(y))(/3(v(x)) - /3(v(y))) dx dy > 0. 

1 J Jk 2N 

Of course, the latter integral may well be (plus) infinity. 







10 


GORO AKAGI, GIULIO SCHIMPERNA, AND ANTONIO SEGATTI 


2.5. Asymptotics of (—A) r and principal eigenvalues as r \ 0. Let us start with the following 
lemma on the behavior of (—A) r as r 0, which will play an important role in the sequel: 

Lemma 2.1. For any £ S(R N ), there holds 

(44) (—A ) r (j) —°> </> strongly in H a (R N ) for any a > 0. 

Proof. The Plancherel identity and the definition of (—A)’' by Fourier transform imply 

( 45 ) ||(-a)-0 - 4 2 h » (RN) = ||(i + Ki 2 r /2 (ier - i)^|| 2 L2(Rf )■ 

Consequently, since f> belongs to the Schwartz class S(R N ) of rapidly decreasing functions, we can 
find a positive L 1 (R^ r )-function g such that 

|(1 + |£| 2 )“(|£r - 1) 2 0 2 (O| < 9(0 for all £ e Rf. 

Then, noting that (1 + |£| 2 )“(|£| 2r — 1) 2< ^ 2 (£) r ' N °> 0 for any £ £ Rjf, we obtain (44) via the dominated 
convergence theorem and (45). □ 

It is also important to recall the following (asymptotic) relation between the H r - and the 
L 2 -norm (see [18] and [26]), 

C(r) 

(46) lim —= IMIi 2 ^ for any v £ |J A Tj0 , 

rG(0,l) 

which is clearly related to Lemma 2.1. 

Next, we shall characterize the behavior of the (weak) fractional Laplacian 2l r as the index r 
goes to 0 (cf. Lemma 2.1). Indeed, we can prove that 2l r tends to the identity operator in a suitable 
way. 

Lemma 2.2. Let {r &} C (0,1) be a sequence with \ 0 as k /* +oo. Let {vk} C Hq be a sequence 
such that Vk £ D(f& rk ) for all k £ N. Moreover, let us assume both {vk} and (24 - k Vk} to be uniformly 
bounded in Hq, respectively. Then, denoting by v the weak limit of Vk in Hq (up to a non-relabeled 
subsequence ofk), we have 

(47) 21 r k Vk - > v weakly in Hq. 


Proof. The boundedness of 2t rfc Vk in Hq entails the existence of w £ Hq such that 

21 rk Vk k ^ +co y w weakly in H 0 , 

up to a non-relabeled subsequence of k (of course, we can assume that this subsequence is extracted 
from the subsequence along which Vk weakly converges to v in Hq). To conclude, we have to prove 
that w = v a.e. in fi. We test the weak convergence 21 rk v>k — > w by tp £ C°°(R N ) with support in fi. 
Recalling Lemma 2.1 and (29), we then have 



wip dx 


lim / 2l rfe Ufc(^da; 

k/' + OO J R AT 


= lim 

k y^+oo 



v k %- k ipdx 


which implies w = v almost everywhere in R w . 



vipdx , 


□ 


The optimal constant of the Poincare-type inequality (20) depends on the first eigenvalue 
Ai(r) of the fractional eigenvalue problem 

(48) v £ X r .o-. 21 r v = Xv in H 0 , 


which also does not mean that 2l r w vanishes outside fl (see §2.3). More precisely, the optimal value 
of cp(r) is given as 1/Ai(r), that is, 


(49) 


INI 


2 

H 0 


< 


1 

Ai(r) 


1^0 


for all v £ X r ,o- 


The first eigenvalue Ai(r) is characterized in the next lemma, which comprises results from [33] and [34] 
and clarifies the spectral properties of the fractional Laplacian. 
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Lemma 2.3 ([33], [34]). Let r G (0,1). Then the fractional eigenvalue problem (48) admits a first 
eigenvalue Ai(r) which is strictly positive, simple, isolated, and can be characterized as 


(50) 


Ai(r) = min 


C(r), 


\ u \h-(R^) = mm 


C(r) NflT(RM) 


“eAf r , 0 ,i|ti||H 0 =i 2 ' ' «£<*•,., 0 \{o} 2 11 11 

Moreover, there exists a unique positive first eigenfunction e\ G X r $ which satisfies || ei || = 1, and 

attains the minimum in (50). 

Finally, denoting with Ai the first eigenvalue of the (standard) Dirichlet problem 

(51) — Av = A v in £1, v = 0 on dfl, 
then for any r G (0,1) there holds 

(52) Ax(r) < AJ. 


For our purposes it will also be extremely important to understand the asymptotics of Ai(r) 
with respect to r \ 0. The following proposition clarifies the situation. In particular, we prove that, 
as r \ 0, the corresponding sequence of first eigenvalues of 2l r tends to 1. Note that, being 21 r an 
approximation of the identity as r \ 0, then the result is not completely unexpected. However, at 
least to our knowledge, this was not observed before. Moreover, due to some lack in compactness 
in the sequence of the first eigenfunctions of 2l r the proof is not straightforward and requires some 
precise estimates. 


Proposition 2.2. Let {r k } C (0,1) satisfy rk \ 0 as k /* +oo. Then we have 

(53) lim Ai(rjfe) = 1. 

+ OO 


Proof. To simplify the notation, for any k G N we denote by A k and Vk the first eigenvalue and 
the (normalized in H 0 ) first eigenfunction of (48), respectively. Then, the couple (A k,Vk) solves the 
eigenvalue problem 

(54) 21 Tk Vk = A kVk in H 0 . 

As a consequence, the only uniform (in k ) estimate on which we can rely is ||'Cfc || /,2 (ir w ) = 1- Hence, 
the only convergence we expect for Vk is a weak L 2 - convergence. Indeed, this is enough to pass to the 
limit in (54) (recall Lemma 2.2). On the other hand, since a priori we cannot exclude that the weak 
limit of Vk is zero, we are not able, at this level, to conclude that A k —> 1. For this reason, we have 
to use a different approach. First of all, thanks to (52) of Lemma 2.3, we only need to prove a lower 
bound for the limit in (53), namely 

(55) lim inf A k > 1. 

k /*+oo 

To this aim, we test (54) by Vk- Using that ||^fc||// 0 = 1 and recalling (29) and (31), we get 

(56) A fc = X k \\vk\\ 2 H 0 = {%r k V k ,Vk)x rk ' 0 

= [ \(-Ay^v k \ 2 dx= [ |£| 2 r ‘MOI 2 d£. 

J R N Jr*) 

For each k G N, set /*,(£) = |Dfc(C)| 2 - Then f k > 0 a.e. in R^ and, by Plancherel’s identity, fk G L 1 (M^ r ) 
with f RN fk(0 d£ — 1 for any k G N. Moreover, by (56) and (52), the sequence of the 2rfc-moments of 
fk is uniformly bounded with respect to k. Namely, we have 


(57) / \e rk fk(Odf = Xk<(Xi) rk - 

J if 

Now, combining this information with the fact that Vk is zero outside H, we can easily deduce that fk 
lies in L°°(R^). Indeed, using |e _ “'^| < 1 together with the Holder inequality, we obtain 


(58) 


fk(0 


< 


( 2 ti -) N 

ini 

( 271 )^ 



e~ ix <v k (x) dx 




\m 

(2n) N 


2 


1 

( 271 -)^ 




Vk(x) dir 


2 


for all ^ G Rf. 
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At this point, we need to estimate the L-norm of f k in a quantitative way. To this purpose, we use 
the following interpolation lemma, whose proof is based on a technique widely used in the kinetic 
theory community (alternatively, one could refer to [38, Lemma 2.1], where a similar inequality is 
shown in a different way): 

Lemma 2.4. Given f £ L 1 (M JV ) (~l L°°(R N ) with f RN |x|“|/(x)| da; < +oo, where a > 0, there holds 

N N 

(59) ii/iili(r") < «Wa)||kr/||^)ii/ii i ^ )I 

where 

k(N, a) = a “+ N (a + N)N~ N + a d N + a 
and d = d(N) is the volume of the unit ball in R w . 

Proof. To simplify the presentation, we set L := ||/||l°°(riv) and M := f RN |x|“|/(x)| dx. We compute, 
for R > 0 to be chosen later, 

(60) [ \f\dx= [ |/| da; + [ \f\dx < dR N L + R~ a M. 

JR N J\x\<R J\x\>R 

Now, optimizing with respect to R the function g(R) := dR N L + R~ a M , we find 

- / otM \ 1 /(-^+ Q! ) 

(61) II/I|li(rN) < mf g(R) = g(R), with R = 

Hence, computing g(R ), we readily get (59). □ 


Continuation of Proof of Proposition 2.2. Applying the lemma, with/ = f k , a = 2r k , |||£| 2rfc /fc||Li(RN) = 
A fc and ||/fe||i~(RJV) < we then get 

(62) IIMlLHRf) < <N, 2r k )\f^ (j^n) > 

which, together with the fact that ||/fc ||z,i(K^ r ) = 1 by construction, implies 


N+2r k ( (2n) N 

(63) Afc > n(N,2r k ) ir ~ k ' — 

Hence, letting k +oo and noting that linifcyi+oo k(N, 2r k ) 
thesis. 



= 1, we get (55), which implies the 

□ 


3. Assumptions and statement of main results 

Let us start with listing our hypotheses on the data of the problem. First of all, we make 
precise the choice of the confining potential by choosing 

(64) W(v) = -\v\ p -\v 2 , p £ (l,oo) \ {2}. 

P 2 

We also set 0(v) := ^|u| p and /3(v) := p{v) = |u| p_1 signu, for brevity. In general, we note the case 
p > 2 as the coercive case, while for p £ (1,2) we speak of a non-coercive potential. Indeed, W is 
unbounded from below in the latter situation. As noted in Introduction, the coercive case includes, 
up to an additive constant, the standard double well potential 

(65) W(v) = i(u 2 - l) 2 . 

Next, for v £ Hq, we introduce the energy functional 

( 66 ) E a (v):=h\v\\ 2 Xit ' 0 + f W{v)dx, 

J Q 

whenever it makes sense. In particular, in the case p > 2, the domain of W, i.e., the set where it takes 
finite values, is given by 

(67) £„ = {u £ Hq such that (u) < +oo} = X„_o fl L P (R N ). 
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Of course, £ a is a Banach space with the natural norm 

(68) IMk ■■= \\v\\ Xai0 + IH| L p (R iY). 

If p < n 2N 2 , then by Sobolev’s embeddings £„ coincides in fact with X a ^ 0 (hence we can use the 
norm of X a $ instead of (68)). On the other hand, in the non-coercive case p £ (1, 2), the functional 
£c may be unbounded from below. In that situation, existence and uniqueness of (global in time) 
solutions still hold. However, when dealing with the singular limits, some difficulties arise from the 
lack of coercivity. 

In what follows, we will look for solutions taking values in the energy space £ a . Correspond¬ 
ingly, we ask that the same regularity is satisfied by the initial datum: 

(69) u 0 £ £ a . 

Next, we introduce a weak (energy) formulation of the system (l)-(2). Here and henceforth, the 
notation C^QO, T]; X) will represent the class of continuous functions on [0,T] in the weak topology 
of a normed space X. 

Definition 3.1. We say that (u,w) is a weak solution to the Cauchy-Dirichlet problem (l)-(4) for 
the fractional Cahn-Hilliard system if, for all T > 0, we have 

(70) u£C w ([0,T};£ a )nC([0,T]-H 0 )nW 1 ’ 2 (0,T;X^ 0 ), 

(71) w£L 2 (Q,T-X sfi )- 

moreover, the couple ( u,w ) satisfies, a.e. in (0 ,T), the following weak formulation of (l)-(2): 

(72) d t u + 2l s ui = 0 in X^ 0 , 

(73) w = + B(u) — u in £' a , 

where B(u ) denotes the bounded linear functional on L P (M. N ) defined by 

(H(u),u) iP(R Ar ) = [ (3(u(x))v{x) dx for v £ L P (R N ), 

J n 

and, finally, the initial condition (3) holds in the following sense: 

(74) u(t) —> uo strongly in Hq and weakly in £„ as t\ 0. 

Correspondingly, we can prove the following existence and uniqueness result: 

Theorem 1. Let s,a £ (0,1), p £ (l,oo) \ {2}, and assume (64) and (69). Then, the fractional 
Cahn-Hilliard system (l)-(4) admits a unique weak solution ( u,w ) in the sense of Definition 3.1, 
which additionally satisfies 

(75) \\/3{u{-,t))\\ 2 Ho < 2 (|Hf)||# 0 + \\u(t)\\ 2 Ho ) for a.e. t £ (0 ,T), 

(76) p(u)£L 2 (0,T-H o ). 

Moreover, u(t) := u(-,t) and¥. a (u(t)) are right-continuous on [0,T) with respect to the strong topolo¬ 
gies of £ a andR, respectively, and the following energy inequality holds true: 

(77) lk(t)||^ 3|0 + <j(u(t)) < 0 for a.e. t £ (0 ,T). 

In particular, if a > s, thenu £ C([0,T];f cr ) andM a (u(t)) is absolutely continuous on [0,T]; moreover, 
the inequality (77) can be replaced by an equality, namely we have 

Mt)\\k, 0 + ^^(“W) = 0 f° r a - e - 1 e ( 0 ; T ). 

A proof of the above result will be carried out in Section 4 by means of time-discretization, 
a-priori estimates, and compactness arguments. It is worth noting that, in addition to (70)-(76), one 
could show that weak solutions satisfy parabolic time-smoothing properties. We will analyze this issue 
in a forthcoming paper, where we also plan to consider a more general class of potentials W. 

Now, we come to the behavior of (families of) solutions as a tends to 0. Then, in the coercive 
case we can prove the following 
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Theorem 2 (from Cahn-Hilliard to porous medium). Let p £ (2, oo), s £ (0,1) and let {or} C (0, 1) 
be such that <Jk \j 0 as k +oo. Moreover, let us given a sequence of initial data and n o £ X' s 0 

satisfying 


(78) 


supE^ (n 0 ,fc) < +oo, n 0 ,fc u 0 in X( 0 . 

feeN 


Let ( Uk,Wk ) denote the corresponding sequence of unique weak solutions to (l)-(4), with a = and 
initial datum Then, there exist a (non-relabeled) subsequence of {k} and a pair of limit functions 
(it, w) such that 


(79) 

(80) 
(81) 
(82) 


Uk —> u weakly star in L°°(0,T; Lq(R n )), 

strongly m 7/(0, T; Lg(R*)) n C([0, T]; *',„), 
weakly in W 1,2 (0, T; X( 0 ), 

Wk —> w weakly in L 2 (0, T ; T Sj o)- 


Moreover, 

(83) (3(u) = w £ L 2 (0, T; T s , 0 ). 
Furthermore, u is a (weak) solution to the fractional porous medium equation 

(84) dtu + 2l s f3(u) = 0 in X' s 0 , a.e. m(0,T), 
with the initial condition n|t=o = no, ie., 

(85) u(t) —> uq strongly in X' s 0 as t \ 0. 


Now we address the non-coercive case p £ (1,2). In this situation (which corresponds to the 
fast diffusion range when dealing with equations of the type (84)) we can prove similar results but 
the analysis requires some extra assumption together with a small modification of the system. First 
of all, we need a compatibility condition between p and the fractional order s in (72). To be precise, 
we ask that 


( 86 ) 


p > 


2N 

N + 2s 


2 *. 


Note that the above condition implies that H S (Q) L p (LI) with compact embedding. Hence X St o is 
compactly embedded in L p (R^). Moreover, we need to modify a bit the energy functional. Namely, 
we replace E CT with the following 

(87) Mu) := + ^\H P LP(n) ~ for v £ £ a = X afi , 


where Ai (a) is the first eigenvalue of (48) with r = a. Note that, in view of Ai (a) <T ^'°> 1 (cf. Prop. 2.2), 
one expects that the contribution of Ai (a) in the limit is idle. The reason why we need to replace E CT 
with Eo- lies in the fact that, thanks to Poincare’s inequality (49), one has 

(88) E ct (i;) > ^|M|£ p(rW) for all v £ X afi , 

and independently of the value of o. Of course, modifying the energy through the choice (87) leads 
correspondingly to a modification of the “original” system (72)-(73), which is now replaced by 

(89) d t u + 2i s w = 0 in X' s 0 , 

(90) w = 2l a u + B(u) — \i(a)u inXf 0 , 

still complemented with the initial condition (74). At fixed tr, existence and uniqueness of energy 
solutions for the modified system (89)-(90) follow from the very same argument given for Theorem 1. 
Moreover, the uniform coercivity provided by (88) permits us to prove a counterpart of Theorem 2 
for the fast-diffusion case in the following theorem: 
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Theorem 3 (from Cahn-Hilliard to fast-diffusion). Let s G (0,1), p G (2*, 2), 2* being given by (86), 
and let {er be as before. Moreover, let us given a sequence of initial data {rto.fc} and uq G X' s 0 
satisfying 

(91) supE CTfc (Mo,fc) <+oo, u 0 ,k u 0 in X' s 0 . 

ke N 

Let ( Uk,Wk ) denote the corresponding sequence of unique weak solutions of (l)-(4), with a = 
W'(v) = fi(v) — \\(<7k)v and initial datum Mo.fc- Then, there exist a (non-relabeled) subsequence of 
{fc} and a pair of limit functions ( u,w ) satisfying (79)-(83). Moreover, u is a (weak) solution to the 
fractional fast-diffusion equation 

(92) d t u + 2l s /3(zi) = 0 in X( 0 , a.e. in (0,T), 
with the initial condition (85). 

Finally, we investigate the behavior of weak solutions to (l)-(4) when the order s of the 
fractional Laplacian in (72) is let tend to 0. In this case, we can prove the following 


Theorem 4 (from Cahn-Hilliard to Allen-Cahn). Let p G (1, oo) \ {2}, o G (0,1) and let {s*,} C (0,1) 
be such that Sk 0 as k +oo. Let also Mo,fe G X a p and u o G Hq satisfy 

(93) sup (E ff (uo,fc) + ||mo, fc|| x 1 ) < 00 > uo,k ^ uq strongly in H 0 . 

k \ s k<°/ 

Let us denote as (uk,Wk) the corresponding sequence of unique weak solutions to (l)-(4) with s = Sk 
and initial datum ito.fc- Then, there exist a (non-relabeled) subsequence of {k} and a pair of limit 
functions ( u, ui) such that 


Uk —t U 


Wk —^ VO 


21 Sk Wk -t w 


weakly star in L°°(0, T\ £ a ), 
strongly in C([0, T]; Hf), 
weakly in W 1 ’ 2 (0, T; £' a ), 
weakly in L 2 (0,T; Hq), 
weakly in T 2 (0,T;£’^). 


Moreover, 

u G C w {[0,T}-,£„) w G L 2 (0,T;ff o ), 

w = 21 a u + B(u) — u, and u is a (weak) solution to the fractional Allen-Cahn equation 
(94) d t u + 2lo.ii + B{u) — u = 0 in E' a , a.e. in (0,T), 

with the initial condition u\ t —o = Uq, i.e., 

u(t) —> uq strongly in Hq and weakly in £ a as t\j 0. 


Theorems 2, 3 and 4 will be proved in Section 5 below. Actually, while the proof of Theorem 4 is 
almost straightforward, taking the limit er 0 will require a careful analysis of the behavior of the 
eigenvalues of 210^ as er*, \ 0. The proof of Theorem 3 is simplified by the fact that we will not start 
from the original system (72)-(73), but rather from its modified version (89)-(90), whose energy is 
bounded from below, uniformly w.r.t. er, thanks to (88). As will be clarified from the proof, this fact 
makes the analysis simpler. Indeed, we will not need to use the results in Prop. 2.2 regarding the 
asymptotics of the first eigenvalue of the fractional Laplacian, which, instead, are essential for proving 
Theorem 2. Actually, we expect that the convergence to the fast diffusion equation could hold, still 
for p G (2*, 2), also for the non-modified fractional Cahn-Hilliard (72)-(73). However, for the moment, 
this remains as an open question. 


4. Existence and uniqueness 

This section is devoted to giving a proof of Theorem 1, i.e., existence and uniqueness of weak 
solutions along with some regularity properties and energy inequalities. In the existence part, the 
basic strategy of a proof is more or less standard, and it consists of time-discretization, a priori 
estimates, compactness arguments to obtain convergence, and Minty’s trick for the identification of 
the limit. On the other hand, we will face some difficulty in deriving the energy inequality (77), 
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particularly, the differentiability of the energy t i—► E CT (it(f)), and for proving the right-continuity of 
u(t) in the strong topology of £ a . These difficulties are due to the simultaneous presence of two 
fractional Laplacians of (possibly) different order. However, the energy inequality seems to be one 
of the minimum requirements for the analysis of singular limits in latter sections and also for the 
investigation of the long-time behavior of global (in time) solutions. Moreover, the right-continuity 
of solutions plays a crucial role to prove the convergence of each global solution as t goes to infinity. 
More details on the long-time behavior of solutions will be discussed in a forthcoming paper. Here, 
it is also worth emphasizing that, in contrast with usual studies on the Cahn-Hilliard equation, our 
argument does not rely on the coercivity of the double-well potential W. In other words, we can treat 
the cases p > 2 and p G (1, 2), simultaneously. This fact potentially means that our proof could also 
be extended to fractional Cahn-Hilliard equations with more general classes of potentials. 


4.1. Uniqueness of solution. We first note that 2l s : X s ^ —> X' s 0 is invertible, i.e., there exists the 
inverse mapping 21J 1 : X' s0 —> X s0 of 2l s . Indeed, the functional on X b> q defined by 

(95) v ^ 0 

is smooth, convex and coercive in X s ^. Hence its Frechet derivative 2l s is invertible. Furthermore, 
one observes by the definition of norms in dual spaces that 

(21 a u, u) Xs 0 = \\u\\ 2 Xs Q = ||2l s u||*, o for all u G X St0 . 

Hence 2l s is a duality mapping between X s q and X' s 0 . 

Now, let (u\,wi) and (u 2 ,w 2 ) be weak solutions of the fractional Cahn-Hilliard equation for 
the same initial data uq. Let u := u± — u 2 and w := W\ — w 2 . Then we have 

(96) dtu + 2l s w) = 0 in Xg 0 , w = + B(ui) — B(u 2 ) — u in £' a . 

Applying 21J 1 : X' s0 —>• X St0 to both sides of the first equation of (96), we have 

d 


dt 


2L 1 u + w = 0 in T 


s,0- 


Test it by u G X s , 0 ^ X' 0 (see (24)) to get 

(97) 

Here we used the fact that 

(98) 


+ (u,w)x a ,o = 0. 


d 

dt' 


u,— a -1 


X,. 


u ) = ( 21, o 2L 1 u, —2L x i 


’ dt 


x a . 


= (t- ll a . _1 «" 2 

2d t 11 s 


1 d 


*•.0 2 dt ]]U 


\% 


since 21,. 1 is a duality mapping between X’ s 0 and X s 0 . Test now the second equation of (96) by u. 
Then, by the monotonicity of /3, it follows that 


(99) 


{w, u) £ „ = \\u\\ 2 x ^ + f W(ui)~ /3{u 2 )) udx- ||u||| 2 (n) 

J 


>INlL,o-Nli3 ( n)- 

We also note that, for any functions u G £„ X' s 0 by (24)) and w G X s> q (°A ££), 

(u,w)x a0 = / u(x)w{x) dx = {vu,u)t 
J RN 


( 100 ) 




Combining (97) and (99) with (100), we deduce that 


1 d 11 ~ 1 1 2 11 ^ 11 2 ^ || ^ || 2 ^ 1 || '-m 2 , || ~ 11 2 


2dt -W^, 0 + INlL.o < ll*lli» ( n) < 0 + c Mx’,o 

for some constant C > 0. Here we used Eluding’s lemma (see [36, Lemma 8]), i.e., for each £ > 0, one 
can take a constant C e > 0 such that 


(101) 


IMU 2 (n) = IMIffo < 4 v \\x„,o + c e \\v\\x> for all v G X a , 0 . 
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Thus we get 


1 d l '- l l2 , 1 || - 1,2 


2dt INI ^o+2 ll “ ll ^^ CNI ^o- 

which along with the fact that n(0) = 0 in X' s 0 implies u(t) = 0 in X' s 0 for all t > 0. Consequently, 
since ui(t) and U 2 {t) belong to H 0 , one can assure that u\(t) = U 2 (t) a.e. in R w for all t £ [0,T], that 
is, u i = U 2 as desired. 

4.2. Time-discretization. Let N £ N and let r = T/N be a time step. In order to construct weak 
solutions, we first carry out the following discretization of (72)-(73): 


( 102 ) 

(103) 


Un U n —1 


+ 2l s u; n = 0 


in X' afi , 


— 21 aUn “b B(Un) Hn— 1 ill £(j 


for n = 1, 2,..., N. In order to show existence of (u n , w n ) satisfying (102)-(103), let us introduce the 
functional F n : £„ —> R given by 

2 


F n(u) := - 


U - U n -1 


r 


C(a) 


X' 


R 2N 


1 11 ' 3 - y\*+^ dxdy + L ^ ^ ~ L Un ~ lU ^ 


for u £ £ ai where 0 denotes the primitive function of 0 such that /3(0) = 0. Then F n is strictly convex, 
coercive and of class C 1 in £„. Therefore one can take a unique minimizer u n £ £ a of F n . Let us 
recall that 21J 1 : X' s 0 —> X St o is the inverse duality mapping, whence Slj 1 coincides with the Frechet 
derivative of the functional 

K(v) :=^\\ v \\%, o for v£X' sfi . 

Correspondingly, K(-) is of class C 1 in X' s 0 ; moreover, 

(104) K\v) = 21” 1 v for v £ X' s 0 . 

Indeed, for each u £ £ a , one has 

K(u) — K{v) < (u — v, ^£i~ 1 u') x u — v) £ for all v £ £ a , 

which implies (104). Thus we have 

7 1 f Un + 2t<7«n + B(u n ) - «n-l = 0 in £' (J . 


(105) 

Setting 


a: 


w n ■— 2L 


1 ( U n U n —\ 


£ X. 


s,0 


and applying to both sides, we obtain (102) along with (103). 

4.3. A priori estimates. Test (105) by u n £ £ a X' s 0 ). Then as in (98) it follows that 
ll u n|| X' 0 — \\ u n— 1 1| x‘ a 0 + 2t (|KII*„, 0 ll Mjl ll LP(Q)) — 2T ( u m u n-l) L 2 (fi) 


< T 


(ll u n||z, 2 (Q) + ||Un-l||ia(f2)) • 


Summing up, we deduce that 

n n 

ll M ™ll*X 0 + 2 X] r (11^11^.0 + ll M .illz,P(fi)) - ll«o||^ i0 + 2 y^ T ll M jllL 2 l'OV 

i =1 3=0 

Recall (101) again to obtain 

n n 

\\ u nfx^ 0 + J2 r (\\ U j\\x„o +2 W U j\\ P LP(Q)) ^ C IK||i2 (n) +<?5>|M 


2 

XL, 
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Hence, due to the discrete Gronwall inequality, we obtain 

N 

(106) 


max || U n \\ x' s 0 T 'y 1 T (ll u n|| Xv'O H Un HlJ>(n)) 


< C. 


n —0 


Test (102) by w n £ X St0 and (103) by (u n — u n _i)/r £ £ a and employ (100) to get 


(107) IKIli 


ll' u «llar < Ti0 ll M n-illEo , ll^-i II EfO ll M ™lli 2 (n) H^n-illEn) 


2 r 


pT 


It 


By summing up, we have 

n | i 111 

(108) E^ll^llko + dKHi.o + -IKIIW) < dMk , 0 + -ho\\ P L* m + oKlli>(n)- 


. 7=1 


Using (101) and (106), assumption (69), and the the boundedness of 2l s : X s $ -X X' s 0 along with 
(102), we then deduce 


N N 

(109) E T \\ W n\\x a , 0 + maX (iKIlL.o + IKII £P(n)) ^ E 


^n—1 


< c. 


Let u T and w T be the piecewise constant interpolants of {u n } and {u> n }, respectively, and let u T be 
the piecewise linear interpolant of {u n }. More precisely, we define u T , u T by 

U-r(t) = U T (t) — " U n X ~ 'U n — i for t £ [t n _i,t n ), 


in X' 


and w T analogously. Then they satisfy 

(110) d t u T + 2l s u) r = 0 

(111) w T = 21 a u T + B(u t ) — u T {- — t) in £' a . 
Moreover, the previous estimates (109) can be rewritten in the form 


(112) f \\w T {t)\\ 2 x dt + SU P (iKWIlLo + IK(*)||£p (R iV)) < C, ( \\d t u T (t)\\ 2 x , dt<C, 
Jo ’ tefo.Ti v ’ v v Jo ’ 


) o te[o,T] 

whence we easily get also 

(H3) 


te[o,T] 


sup I \\u T (t)\\% o 0 + ||u/^ l|P 


LP 


< C. 


Note that 

{B{u T {t)),(j)) LP(R N ) = J P(u T (x,t))<l)(x) dx < ||ur(i)ll^( 1 f 2 )II^IU' , (» JV ) for a11 0 e L P (R N ), 
which implies 


(114) 


sup ||B(u T (t))|| ip / (RW) < sup \\u T (t)\\ P LP(n) <C. 
te[o,T] te[o,T] 


4.4. Convergence as t — > 0. From the estimates established so far, one can take a (non-relabeled) 
subsequence of r —> 0 (equivalently, N —> oo) such that 


(115) 

w T —^ w 

weakly in L 2 (0, T ; X s $), 

(116) 

21 s w T -X 21 s w 

weakly in L 2 (0, T; X' afi ), 

(117) 

u T —^ u 

weakly star in L°°(0, T ; £„), 

(118) 

%i a u T -X 21 a u 

weakly in L 2 (0, T; X' afi ), 

(119) 

B{u t (•)) -x x 

weakly star in L°°(0,T;L p (M. N )) 

(120) 

u T —^ u 

weakly star in L°°(0, T ; £„), 

(121) 

d t u T -x d t u 

weakly in T 2 (0, T; X' s q). 
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Combining these facts with the compact embeddings £ a H 0 s- X' s 0 , and using the Aubin-Lions- 
Simon compactness lemma (see [36, Theorem 5]), one can verify that 

(122) u T —> u strongly in C([0, T]; H 0 ). 

Then u belongs to C([0, T]; H 0 ) fl C w ([ 0, T]; £ a ) as well. Moreover, we observe by (109) that 


(123) \\u T {t) - u T (t)\\x' 0 = u„ 


t t"n— 1 . in i 

U n + U n —i 


t n — t (109) 

= — -||u n - u n -i\\x', < Cy/r for all t £ [t n -i, t n ), 

T s,0 

which along with (122) yields 

sup \\u T {t)-u{t)\\ x > < sup \\u T (t) - u T (t)\\ x > + sup \\u T {t)-u{t)\\ x > ->0. 
te[o,T] ’ te[o,T] ' te[o,T] 

Moreover, exploiting (101) and (112), for any e > 0, one can take C e > 0 such that 

\\u T {t) - u(t)\\ Ho < e\\u T (t) - u{t)\\ x<r 0 + C e \\u T (t) - u(t)\\ x > o 

<eC + C E sup ||u T (r) — u(r)\\x’ for any t£[0,T], 

re[0,T] 

whence there follows that 

(124) u T —> u strongly in L°°(0,T; H q ), 

(125) u T (t) —> u(t) strongly in Hq for all t € [0,T], 

(126) weakly in £ G for all t£[0,T]. 

Furthermore, noting that 

IK(i-r) - u{t)\\ X ’ s o 

< || U r (t-T) -U r {t-T) \\x’+ II U T {t-T) - U T (t)\\ X > + || U T (t) ~ u(t) || X' 


< 


C\fr + f \\d r u T (r)\\x' 0 &r+\\u T {t) 

J t — T 


Wxt 


for all t £ [t, T], we also deduce from (112) and (122) that 

sup \\u T (t - t) - u(t)\\x' < Cy/r + sup \\u T {t)-u{t)\\x' n -t 0 as r -A 0. 
te[r,T] ’ 4G[0,T] 

As in (124), one can further obtain 

(127) u T (-— r) —>u strongly in L°°(0,T; Hq), 

(128) u T (t — t) — > u{t) strongly in H 0 for all t > 0. 

In particular, (128) also implies w = 2l CT u + x — u in £' a . 

We next verify that x = B(u) by using Minty’s trick. To this end, we observe that 

{B(u T (t)),U T (t)) L pQ l N ) = (w T (t),U T (t)) - \\u r {t)\\ 2 x at0 + (u T (t - t ), U T (t)). 

Hence, using (115), (124) and (128), we obtain 

limsup / (B(u T (t)),u T (t)) LP mN\ dt 

T-> 0 JO 


(129) 


= lim [ (w T (t),u T (t)) dt — liminf [ |K(f)||^ d t 

T ^°J o T ^° Jo 

+ lim / (u T (t — t), u t (t)) dt 

<[ (w{t),u(t)) dt - I \\u(t)\\ 2 x a0 dt+ ( (u(t),u(t)) dt 
Jo Jo ’ Jo 

= / (x(*).w(*))lj>(r") dt, 

Jo 
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which along with (117), (119) and the maximal monotonicity of the operator u H»• B(u) from L P (M. N ) 
to L p (R w ) yields x = B(u) in L p (1^). Thus (u, w) is a weak solution of (l)-(4). 


4.5. Energy inequalities. In order to derive (75), let us fix t £ (0, T ) at which (73) holds and define 
the Yosida approximation f} E : R —> R of /3, i.e., (3 e (r) := (i— j E (r))/e = /3(j E (r )) for r £ R, where j e 
stands for the resolvent of /3 defined by j £ (r) := (1 + e / 9) _1 (r). Then, since /3 e is Lipschitz continuous 
and /3 e (0) = 0, one can observe that (3 E (u(-)) £ £ a if u £ £ a - Hence, we can test (73) by /3 e (u(-,t)) to 
get 

(21 a u(t), /3 e (u(-, t))) Xa o + Pe(u(-, t))) £a = (w(t) + u(t), /3 e (m(-, t))). 

Here we note as in (42) that 

(2l a u(t),p e (u(-,t)))x„' 0 > 0. 

Moreover, by the definition of Yosida approximation and the monotonicity of /?, we infer that 
(B(u(t)),f3 E (u(-,t))) £ir = J /3(u(x,t))p s (u(x,t)) dx > \\/3 e (u(-,t))\\ 2 L 2 {n) = \\Pe(u{-,t))\\ 2 Ho . 

Thus we obtain 


(130) \\P e (u(-,t))\\ Ho < |H<)||ff 0 + ||u(i)||u 0 for a.e. t £ (0,T), 

which implies 

Pe(u(-, t)) ~^b t weakly in H 0 as e —> 0 

for some b t £ H 0 - On the other hand, let us recall that fl e (r) = f3(j e (r)). Moreover, since j e is 
non-expansive (i.e., Lipschitz continuous with the Lipschitz constant 1) and j E ( 0) = 0, one can easily 
check that 

\\je(u(-,t))\\ £a < \\u(t)\\ s<7 , 

which yields 

j e (u(-,t)) —> u(t) weakly in £ a and strongly in Hq, 
as e —> 0. Here we also used the fact that 


\j E (u(x,t)) — u(x,t)\ < e\/3 E (u(x,t))\ < e\f3(u(x,t))\ for a.e. x £ H. 

Therefore, by virtue of the demiclosedness of the maximal monotone operator u i— > f3(u(-)) in Hq x Hq, 
we conclude that bt = /3(u(-,t)) a.e. in H. Moreover, (75) follows from the weak lower-semicontinuity 
of the norm || • ||# 0 in Hq. Furthermore, integrating (75) in time, we obtain j3(u) £ L 2 (0,T; Hq), as u 
and w belong to L 2 (0,T\Hq). 

We shall finally prove that t i —> E a (u(t)) is differentiable a.e. in (0 ,T) and derive the energy 
inequality, namely 

(131) \\w(t)\\x s , 0 + a{u(t)) < 0 for a.e. t £ (0, T), 

which can be also rewritten as 

(132) (d t u(t),w(t)) Xs 0 > a(u(t)) for a.e. t £ (0,T). 

Moreover, the right-continuity of the function 1 1 —> u(t) in the strong topology of £ a will also follow as 
a by-product of our argument. 

Remark 4.1. Before proceeding with a proof, it is worth stressing that, differently from what happens 
in the non-fractional case, the differentiability of E a(u(t)) and the energy inequality (131) are not 
straightforward. Indeed, the energy functional E CT is smooth in £ a but non-convex. Hence, if one 
attempts to apply a standard chain-rule to E CT and u(t), the differentiability of u(t) in the strong 
topology of £ a is needed. However, t u(t) turns out to be differentiable only in the weaker space 
X' s q. When dealing with the standard Cahn-Hilliard equation (i.e., for s = a = 1), this problem 
may be overcome by rewriting the energy functional corresponding to E CT as the sum of a convex 
and of a concave part and by applying a generalized chain-rule for convex but (possibly) non-smooth 
functionals (see, e.g., [8]). However, in the present case, this kind of procedure seems to work only 
when a > s. We shall give the highlights of a proof of this fact in Subsec. 4.6 below. 
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In order to show (131), we start with noting that, from (108) and interpolation, there follows 
J Q \\wr{r)\\% a 0 dr + Ca(u T (t)) - ^||ur(i)||| 2 (n) < C CT (u 0 ) - ^||uo||| 2( n) for 0 < t < T, 
where C CT (-) denotes the convex functional of class C 1 on £ a given by 

^(n) = ^IMIa^o + ~"IMIz,p(n) for vsE,,. 

z p 

Using the convergence relations obtained so far and the weak lower semicontinuity of C CT (-) in £ a , we 
deduce that 

(133) [ ||w(r)||^ s o d r + E a (u(t)) < E a (u(0)) for all t G [0, T\. 

J o 

From the uniqueness of the solution and the fact that u{t) £ £ a for all t £ [0, T], one can also derive 

(134) ^ ||in(r)||^- s Q dr + E <7 (u(r)) — E a (u(t)) < 0 for all 0 <t<r<T, 

which also implies that E ct (m(-)) is nonincreasing on [0,T], whence it is differentiable a.e. in (0,T). 
Since u £ C([0,T]; H 0 ) fl C,„([0, T\; £ a ) and C CT (-) is weakly lower semicontinuous in £ a , E CT (u(-)) is 
right-continuous on [0,T), i.e., E (T (m(t)) —> E a (u(t)) as t \ t. Then, the same property holds for 
||u(-)||i„ 0 and II m (‘)IIlp(rn) dy the weak lower semicontinuity of the norms. Therefore due to the 
uniform convexity of and L p ( R w ), we can also verify that u is right-continuous on [0,T) in the 
strong topology of £ a . 

Furthermore, let t belong to the set 

X \= [t £ [0, T]: E (T (u(-)) is differentiable at t, and t is a Lebesgue point of ||in 

Then (0,T) \ X has zero Lebesgue measure. Dividing both sides of (134) by r — t > 0 and passing to 
the limit as r \ t, we obtain (131). 

4.6. Energy equality. We prove here that, under the condition a > s, u belongs to C([0,T]; £ a ), 
the energy E (T (M(t)) is absolutely continuous on [0, T], and the inequality (77) can be replaced by the 
following energy identity: 

(135) ||w(i)|& s o + ^E a (u(t)) = 0 for a.e. t £ (0, T). 

The key tool in order to get (135) is the following chain-rule formula, which can be proved by adapting 
the argument given in [30, Lemma 4.1]: 

Lemma 4.1. Let be a Hilbert triple and let T : H —> (—oo, +oo] be a convex, proper and 

lower semicontinuous functional. Moreover, let us assume that, for some fci > 0, k 2 > 0, there holds 

(136) ^(u) > fci||n||^ — fc 2 

Denote with A the subdifferential o/4/ with respect to the scalar product ofH, and consider, forT > 0, 
v £ W 1,2 ( 0, T ; V') D L 2 (0, T ; V) and p £ L 2 (0,T- V) with rj(t) £ Av(t) for a.e. t £ (0, T). Then, the 
function 1 1 —► 4 >{y{t)) is absolutely continuous in [0, T]. Moreover, 

(137) f (dtv(r ), t)(t)}v dr = 4/(n(t)) — 4/(n(r)) for all 0 < r < t < T. 

J r 

We apply the above Lemma with the following choices: 

. H = Ho,V = X sfi , see (24). 

• = ±\\v \\ 2 x<70 + J n fj{v(x))dx. 

Then, clearly, is proper, lower semicontinuous and convex. Moreover, thanks to the fractional 
Poincare inequality (20), it satisfies the coercivity assumption (136). Now, let u, w be the solution 
given by Theorem 1. Then, by Definition 3.1, we have 

u £ W 1 ’ 2 ^, T; V') and w £ L 2 ( 0, T; V). 

On the other hand, being <r > s, we also have 

u £ L 2 (0, T ; V). 
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Thus, setting 77 := w + u, it follows that 77 £ L 2 (0,T; V); moreover, thanks to equation (73), 77 (f) £ 
9T(it(t)) for a.e. t £ (0, T). Hence, (135) follows from Lemma 4.1. 

Remark 4.2. When er < s, we do not know whether rj = w+u £ L 2 (0, T; V) = L 2 (0, T; X Si o), because 
X a ,o is not included into X s fi- Hence, the assumptions of Lemma 4.1 are not necessarily satisfied and 
the energy identity (135) remains as an open issue in this case. 

5. Singular limits 


5.1. Uniform estimates. Before proving our results on singular limits as er —> 0 or s —> 0, we 
first establish uniform estimates with respect to s and er for the unique weak solution ( u , w) to the 
fractional Cahn-Hilliard system (l)-(4). To this end, let us recall the inequality (75) of Theorem 1, 

(138) \W(u(-,t))\\ 2 Ho < 2 (IKOII/io + IH*)IIh 0 ) fora.e. t £ (0,T). 

Moreover, we also recall the energy inequality (77), 


^ K(u(t )) + \\w(t)\\ 2 Xs 0 < 0 
Then, integrating both sides in t we get 


(139) 


sup 

te[o,T] 






W(u(x,t))dx 


for a.e. t £ (0, T). 


+ IMlL 2 ( 0 ,T;A, ii0 ) < Q0M«o)). 


Here and henceforth, Q(-) denotes a computable nonnegative-valued function which is monotonely 
increasing in its argument(s) and may vary from line to line. In particular, the expression of Q may 
depend on p, T and |f2|; however, it is always independent both of a and of s. 

In case p > 2: since W is coercive, it follows immediately from (139) that 

(140) Nil + IMIz,=o(0,T;LP(fi)) + ll u 'lli 2 (0,T;A’ a|O ) — Q(Ecr(uo))) 

which along with the boundedness of 2l s : X St0 —> X' 0 and a comparison of terms in (72) gives 

(141) ll^t u lli 2 (0,T;AT' 0 ) = ||2l s «l||| ( a(o ) T;^' 0 ) ^ Q(E(t(uo))- 
Moreover, by (114), 

(142) II^( m )IIt“( 0 ,T;Lp'(R , v ‘)) — ll M llL°°( 0 ,T;LP(n)) — Q(E<r( u o))- 
On the other hand, combining (140) with (138), we get 

(143) ll^( u )lli 2 (0,T;ILo) — Q(®'cr(uo))- 
In turn, this estimate clearly implies that 

(144) \\W\u)\\1 2{0 ^h 0 )<Q^M). 

Hence by (73), we find that 

(21 <7 u(t),(/)) £<T (7 = ( w(t ) - B(u(t)) +u{t),<f>) £tr 

= / (w(x, t) — /3(u(x, t)) + u(x, f)) 4>[x) dx 

Jn 

< (lk(f)||ff 0 + ||/3(u(-,t))[| ffo + |Hf)||ff 0 ) UWh 0 for all $ £ @n(R N ). 


Since ^(R^) is dense in ff 0 , one has 

(145) ll^<T U llz,2(0,T;ff o ): < Q(E CT (7lo))) 

where 2l CT : Hq — > Hq stands for the ifo-fractional Laplacian with domain U(2t cr ) C Hq (see §2.4) and 
2l CT u(f) : H 0 —> R is the unique bounded linear extension onto H 0 of the functional %i a u(t) : X a0 —y R. 

In case 1 < p < 2 and a £ (0,1) is fixed: (the argument below is still available for p > 2 as 
well) by applying 21J 1 to both sides of (72) and by utilizing (73), we have 

217 1 ( dtu(t )) + 21 a u(t) + B{u{t)) — u(t) = 0 in £' a , 0 <t <T. 


FRACTIONAL CAHN-HILLIARD, ALLEN-CAHN AND POROUS MEDIUM EQUATIONS 


23 


Test it by u(t) £ £ a . It follows that 

\^\\ u { t )fx' s0 + IHOIlL.o + IW*)IIL> ( n) = IK*)lll 2 (f!) for a.e. o < t < T. 

Set X = X r o or X = £ r := X r o n L p ( R w ) for a fixed constant r £ (s, 1) or X = Ug(f2). 


Then X is continuously embedded in X St o uniformly for s —> 0. More precisely, there exists a constant 
Co > 0 independent of s —> 0 such that 

(146) IMU ..0 < Colkllx for all v € X. 

Indeed, as in [18, Proof of Proposition 2.1], one can verify that, for all v £ X, 


C(s ) 


Hx)-y(y) I s 

' R 2N \x — y\ N + 2s 


dx d y 


C(s) 


l / 

Jr n U. 


v(x)-v(y) | 2 

aT dx dy 


C(s ) 


{iel N : \x-y\>l} \ x y\ N+ 

\v{x)-v(y)\ 2 


E» j{xeR N : \x-y\<l} \ x ~ y\ N+2s 


dx dy 


< 


C( s ) 10JV-I111 ,||2 


-is 7 ' 


lli 2 (n) 


C(r] \\ V \\Xr,o 
s^- 1 ! C(s) 


for X = or £, 


-r-^IIVwll^n) fOT X = H^Sl). 

Here stands for the surface area of a unit sphere in R w . Finally, exploit the asymptotics (12) 

of C(r) as r \ 0 to obtain (146). Moreover, (146) yields 

(147) IICIIa-' < Co|KIU; i0 for all C e 

which particularly gives 

(148) IMU' < C 0 ||u||^/ Q for all u £ H 0 (~ H' 0 ). 

Indeed, for any (i£lc X s q, one Hnds that 

(CA)x = «,<t>)x..o < IICIU', 0 |Mk, 0 < CollClUiJI^IU for C £ < 0 , 

which gives (147). On the other hand, from the dense and compact embeddings X a> o Ho (— H' 0 ) ^ 
X' along with Ehrling’s compactness lemma [36, Lemma 8], for any e > 0 there exists a positive 
constant C £jfT , which is independent of s but may depend on cr, such that 

(149) IMIi 2 ( n) = Hir 0 < e|MI*„, 0 + C eA v tx' fOT all v £ X afi . 

Therefore we deduce that 

^IK*)ll*,' i0 + l\Wm 2 X<r ,o + IK*)lliP(fi) < C^vM^Wx' f01 ' a.e. 0 < t < T. 

The integration of both sides over (0 ,t) along with (148) yields 


ll«(*)ll 


X' 


(hMIlL.o + IH r )lli,p(n)) dr ^ c <r (|M |^ 0 + f ||w(t)||x' dr 


for some constant C a > 0 (depending on cr). Hence, exploiting Gronwall’s inequality, we obtain 


(150) 


sup im 

te[o,T] 




(lKi)HL,o + IKf)llk (n) ) dt < Q(c*, \\M%J- 


Apply (149) to (139) and employ (150). Then we obtain (140) with a bound depending on C a , ||uo|| x' B 
and E CT (uo). Furthermore, relations analogous to (141)-(145) also follow with similar bounds. More 
precisely, one deduces that 

(151) IMll=°(0,T;A’„ |O ) + IMIz,“(0,T;LP(fi)) + ll/4( u )lll 2 (0,T;ffo) + II W (u) ||| 2 (o i T;Ef 0 ) 

+ l|21cr u l|_L=o(o,T;A^ „) — Q(dlcri ||rto||^ s ' 0 ,E cr (rto))- 

Moreover, we also have 

ll^t' u lli 2 (0,T;A'^ 0 ) + l|24 u ’lli 2 (0,T;.y; 0 ) + ll' u, lli 2 (0.T;A s ,o) — ll'W'O || AC' 0 J E CT (tto)). 
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Hence by virtue of (147) and Poincare’s inequality (49) along with Proposition 2.2, it follows that 

(152) \\d t u\\l^ T , xl) + |MI Z, 2 (0,T;A') + \\ w \\l 2 (0,T-Ho) <Q(C a ,\\u 0 f Ko ,E a (u 0 )). 

Finally, by (114), it holds that 

(153) l|B(«)ll i - (0l r ■LP 1 (R w )) < Q{C„, ||wo||^ 0 ,E<t(uo)). 

5.2. Limit of fractional Laplacian in Bochner spaces. In this section, we shall generalize Lemma 
2.2 for later use of proving the convergence of 21 Tk Uk as \ 0 in an appropriate Bochner space. 
Throughout this subsection, we use the notation Xg and Xgx, even for /3 > 1 in the following sense 

Xg := HP(R n ) = {«£ S(R n )' : (l + |f | 2 ) 0/2 u(f) £ L 2 (Rf)} , 

Xp t Q := {u £ Xg : u = 0 in WL N \ fl} for ft > 1. 

Then one finds that 

Xg -A -a X 0 = L 2 ^), Xg, o -A Xy, 0 ^ X 0 , 0 = H 0 if p > 7 > 0 
with continuous densely defined canonical injections. Hence we also have dual relations, 

L 2 (R N y -A -A X'g, H'o °A xy 0 °A Xg 0 if P > 7 > 0 

densely and continuously. For each u £ Xg,o and / 3 ,y > 0, one can define T(u) £ X’’ by 

(T(u), <j>)x y := f u(x)(/)(x)dx for <j> £ X 7 . 

Jn 

Then T : Xgg -a X’ is continuous due to the continuous embeddings described above. Hence Xgx, is 
continuously embedded in X ’’ by T. From now on, we simply write u instead of T(u) if no confusion 
may arise. 


Lemma 5.1. Let u and f be integrable functions of (0, T ) with values in X’g 0 for some constant ft > 0 
satisfying ft Y n — 1/2 with n £ N. Let {r^} be a sequence in (0 , ft) such that \ 0 as k —> oo and 
consider a sequence {uk} of strongly measurable functions in (0, T) with values in X rk ,o. In addition, 
assume that 


(154) 

(155) 


{u k (t),(pk) Xp (fit) d t 


(&r k Uk(t),<p) Xrk ' 0 <l>(t) dt 


J {u(t),<p) Xpi0 ^>(t) dt, 

[ d>(t)dt 

Jo 


for any p £ ^(1^), <p £ Cg°(0,T) and (pk £ Xg satisfying tpk —> P strongly in Xg. Then it holds 
that u(t ) = f(t) in Xg 0 for a.e. t £ (0,T). 


Proof. From the continuous embeddings X rk , 0 A L 2 (R N ) ~ (L 2 (R Ar ))' A Xg, we note that the 
function u k is strongly measurable with values in Xg as well. For any tp £ ®si(R W ) and (j) £ C^°( 0, T), 
we observe by (154) that 

rT 


(156) 


(21 r h u k {t),(p) x <j>{t) dt 


C(r k ) 


2 

t-T 




i(~A) rk ip,u k it)) (fit) dt 


= (Uk{t ), (-A ) rk (p) Xf> (fit) dt (1 -a 4) J (u(t), <p) Xf> 0 (fit) dt. 

Here we used the fact that (— A) rk ip £ Xg and (—A ) rk p -A <p strongly in Xg (uniformly in t) as k — > oo 
by Lemma 2.1. By virtue of (155), the left-hand side of (156) converges as follows: 


(21 T h u k {t),<p) x <f{t) dt 


(155) 


(£&)>?)X/,,0 (fit) dt. 
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Thus we obtain 


J (£(*)> <£(*) dt = J <“(*). 0W di - 


Recall that %(l w ) is dense in dlgg if /? ^ n — 1/2 for n G N (see [25, Theorems 11.4 and 11.1, 
Chap. I]). Hence, from the arbitrariness of p G £>q(R n ), we have 

[ (£(i) — u(t)) <f>(t) dt = 0 in X'p o for all (f> G C£°(0,T). 

Jo 

Finally, applying du Bois-Reymond’ lemma for Bochner integrals, we conclude that £(i) = u(t) in 
X'g 0 for a.e. t G (0,T). □ 


Remark 5.1. (i) All the assumptions of Lemma 5.1 can be proved to hold whenever Uk —t u 

weakly in L l {d,T] Xg) and 21 rk Uk{-) —> £ weakly in L 1 (0,T;A^ 0 ) as k — > oo for some u G 
i 1 (0,T; Xg) and £ G T 1 (0 ,T;Xp 0 ). Indeed, the product of test functions pkj> converges to 
p4> strongly in T°°(0, T; Xg) and p<f> belongs to L°°(0,T; Xg^), since pk —> p strongly in Xg 
with p G @n(R N ) C Xp t Q and (j) G Cg°(0, T). Moreover, since Xg 0 C Xg. one observes that 
Xg H- Xg f) . whence u belongs to L 1 (0,T; Xg 0 ). 

(ii) One can also derive a similar result for sequences independent of t. More precisely, if Uk G 
X rkj o, Uk —> u weakly in Xg. 2 i rk Uk G Xg (j and 21 Tk Uk —> £ weakly in Xg Q , then £ = u in Xg Q . 
Indeed, set Uk(-) = Uk■ Then Uk —> v(-) = u weakly in L p (0,T; Xg) and 21 rk Uk{-) = 21 rk Uk 
converges to rj(-) = £ weakly in L p (0,T\Xg 0 ) for any p G [l,oo); hence, all the assumptions 
of Lemma 5.1 hold true by (i) above. 


5.3. Proof of Theorem 2. Let {cr fe } be a sequence in (0,1) such that Ok \ 0 and let (uk,Wk) be 
the family of weak solutions to 

(157) d t u k + 21 s w k = 0 in X' s0 , 

(158) w k = %a k Uk + B(u k ) - u k in £' ak 


with uk( 0) = uo,k- Then recalling uniform estimates (140)-(145) in §5.1 along with hypothesis (78), 
one can take weak limits u G L°°(0, T; L^(R N )) r\W 1,2 (0, T; X^ 0 ), w G L 2 (0,T; T Sj0 ), j3 G L 2 (0,T; H 0 ) 
and £ G L 2 (0, T; L/g) such that, up to a (non-relabeled) subsequence of {fc}, 


(159) 

Uk —^ U 

weakly star in L°°(0, T; Lq(M. n )), 

(160) 


weakly in IT 1 ’ 2 (0, T ; X' s0 ), 

(161) 

Wk —> w 

weakly in L 2 (0, T; X s ^), 

(162) 

21 s w k -A- 2 t s ui 

weakly in L 2 (Q,T-,X' s0 ), 

(163) 

/3{u k ) —t /3 

weakly in L 2 (Q,T\H 0 ), 

(164) 


weakly in L 2 (Q,T\H q ). 


It follows immediately that d t u + 2l s w = 0 in X' s 0 . Applying Lemma 5.1 with any (3 > 0 to 21 ak Uk{t) 
and Uk{t) along with the weak convergence relations (159) and (164), we obtain £(t) = u(t) in Xg 0 
for a.e. t G (0,T). Moreover, since £(t) and u(t) lie in H 0 , which is dense in Xg 0 , we see that £ = u 
a.e. in Q x (0,T). For all (f> G ^(R^) and p G Cq°(0,T), it follows from (158) that 

[ (w k (t), <p) 4>{t) dt = [ (f2L ak Uk(t),p)<t>(t) dt 

Jo Jo 

+ [ (P(u k (-,t)),<p)<t>{t) dt- [ (u k {t),p)<j>(t) dt. 

Jo Jo 

Passing to the limit as k /* oo, we obtain 


(w(t), p) <j>{t) dt = (£(f), p) (j>(t) dt 


(P(t),p) <Kt) dt- [ (u(t),p)^>(t) dt, 
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which together with the density of in Hq and the arbitrariness of </> £ C{j°(0,T) implies that 

w = f; + 0 — u in Ho, a.e. in (0, T). 

Thus we obtain w[t) = 0{t) in H 0 for a.e. t £ (0, T), or, in other words, 

(165) w = 0 a.e. in 12 x (0,T). 


For each t £ [0, T], since {itfc(t)} is bounded in iJ 0 and i/ 0 is compactly embedded in X' s 0 , the 
sequence {wfc(t)} is precompact in X a0 . Moreover, t u k {t) is equicontinuous on [0,T] with values 
in X' s 0 for k £ N. Therefore, thanks to Ascoli’s lemma, we infer that 


(166) 


u k ^u strongly in C([0, T\\X' s0 ). 


Since uo,k —> uo strongly in X' a 0 by assumption, one can check that u(t) —> uq strongly in X' s 0 as 
t \ 0. In particular, w(0) = uq. 

Now, the major task is to identify the limit 0 as 0 (u), namely proving that 0 = 0{u) a.e. in 
x (0,T). To this end, we shall use Minty’s trick, i.e., we claim that 


(167) lirnsup / / 0{u k )u k decdt < 

fc/^+oo Jo JR N 

Actually, testing (158) by u k , we find that 

[ [ 0{u k )u k dxdt = - [ \\u k (t)\\x„ o dt+ f 

I o Jr n Jo k ' Jo 

Taking the limsup fc * +OG of both sides, we have 
(169) lirnsup [ [ 0{u k )u k dxdt 


lo Jr n 


0u dx d t. 


(168) 


ll u fc(^)ll l 2 (r n ) ' 


w k u k dx dt. 


lo Jr n 


k/'+oo Jo J I 

r 1 

< lirnsup / 
fc/^ + OO J 0 


M feWllL2 (RN) - \\U k {t)\\ Xc 


'fc.o) 


dt 


wu dx dt. 


lo Jwl n 


In particular, in order to take the limit of the last integral, we used (161) together with (166), and 
observed that 


wu dx dt. 


/ o J«. N 


(170) [ ( w k u k dxdt= j (uk(t),w k (t))x a , 0 dt -A f (u(t),w(t)) Xs 0 dt 

Jo Jr n Jo Jo 

The Poincare inequality (49) gives 

(171) D k (t) := - \\uk(t)\\ 2 Xrrk 0 < W u k(t)\\x„ h , 0 for a - e - t S (0, T). 

Thus, recalling Proposition 2.2 and the energy estimate (139), we conclude that 

lirnsup f D k (t)dt<0. 

oo JO 


Then (167) follows from the above along with (165), (169) and the fact that u = 0 outside 12. Therefore 
thanks to (159), (163), and the maximal monotonicity of the mapping u H > 0(u) in L 2 { R w ) x L 2 ( l w ), 
one deduces that 0{u) = 0 a.e. in R w x (0, T). In particular, 0 = (3(u) vanishes outside 12. Hence by 
(165) together with the fact that w = 0 in R N \ 12, one obtains 

(172) /3(it) = 0 = w a.e. in x (0,T), 

which also implies 0(u) =w£ L 2 (0,T; X s $). Consequently, u solves for almost any t £ (0, T) 

d t u + 21 3 0(u) = 0 in X a 0 . 

Note that, as a consequence of the procedure, recalling (159) again, we also get 

(173) u k —y u strongly in L p ( 0, T; Lq(M. n )) 
by utilizing the uniform convexity of L^M^). 
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5.4. Proof of Theorem 3. We first remark that, as in the Riesz representation theorem for standard 
Lebesgue spaces, one can also identify the dual space (Lq(R n ))' of Lq(WL n ) with Lg (1 w ), where 
q € (1, oo) and q' := q/(q — 1). 

Let (7fe \ 0 and let (uk, Wk) be the family of weak solutions to 


(174) d t u k + 2l s w; fc = 0 in X^ 0 , 

(175) Wk = ^ia k Uk + B(uk) - XkUk in X' akfi , 

with Uk( 0) = uo,k, where \ k := Xi(<Jk) denotes the first eigenvalue of (48) with r replaced by o>. As 
in §5.1, (formally) test (174) by w k and (175) by d t uk to get 

IMf)ll*., 0 + < 0 for a.e. t e (o, t), 

where E CTfc : X aki o — > R is defined as in (87). Indeed, the energy inequality above can be rigorously 
derived as in the proof of Theorem 1. Integrate both sides over (0, £). It follows that 

ll' u; fe|lL 2 (0,T;A s ,o) + (ufc(t)) < E ak (uo,k)- 


By (88), we have 

ll' u ’fc|lL 2 (0,T;A s ,o) + H Mfc llz / °°(0,T;Lg(R iv )) — Q(^k i u 0,k)), 

which also implies 


ll^t M fc|ll, 2 (0,T;A' iO ) + ll^s u, fc|li 2 (0,T;^ 0 ) + ll^( u fc)ll^ oo ^ 0 T . L p'( R iV)) — ( 3(®'°'fc ( u 0,k))- 

As in (145), by (175) and estimates above along with the fact that 1 < p < 2 (i.e., p' > 2), we can 
take the unique bounded linear extension ( H< Tk Uk{t) : Lq (R^) —> 


onto Ll 


2fo-fcWfe(£)| i p'^ Rjv ^ : X Gk $ 


n Ll 


such that, by the identification (Lq 


of the functional 


— 


a k Uk(t) 


L p 0 ( R w ) 


< ||xt7fc (*)||z,g(]RJV) + Afc||u fc (£)|| i P( R JV) + \\/3(u k (•,£)) II Lg(RW) 

< C (|K(£)||il 0 + Afc||ujfe(i)|| i g (R i\i' ) + IM^II^V)) 


for some constant C > 0 independent of k and t. We shall simply write 21 ak Uk instead of 2l crfc zifc(-) 
below. Thus we obtain 

|| 2 l(Tfc' U fe|li 2 (g j j’.j / P( R W^ < Q0&a k (uo,k))- 

Therefore, there exist weak limits u, w, j3 and £ such that, up to a (non-relabeled) subsequence, 


(176) 

Uk —l u 

weakly star in L°°(0, T ; Lq(R 7 ' 

(177) 


weakly in W 1,2 (0, T; X' >0 ), 

(178) 

Wk —^ U1 

weakly in L 2 (0, T ; X sfi ), 

(179) 

P{uk) -*■ P 

weakly star in L°°(0 ,T;Lq (R 

(180) 

2l(j fc Uk ^ £ 

weakly in L 2 (0, T ; Lq(M. n )). 


Moreover, apply Ascoli’s compactness lemma along with the compact embedding Lq( R w ) ~ (Lq (R^))' 
X' sG (see (86)) to get 

(181) u k ^u strongly in C([0,T}; X' s0 ). 

To prove £ = u, we use Lemma 5.1. Indeed, choose ft sufhciently large so that Lg(R w ) is densely and 
continuously embedded in Xp ( X' 3 0 ). Then the weak (star) convergence of u k (cf. (176)) and that 

of 21 <j k Uk (cf. (180)) suffice to apply the lemma (see (i) of Remark 5.1) and obtain the conclusion. As 
Afc —> 1 by Prop. 2.2, we then arrive at 

d t u + 21 s w = 0 in X' %0 , w = /3 in Xp 0 . 
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Hence w = j3 in fl x (0, T). It remains to prove that ft = /3(u) in x (0, T). By Poincare’s inequality 
(49) and (153), we note that 

r T r rT p rT r rT r 


. A™ A-t- 


Jo Jr n Jo Jr n Jo Jr n Jo Jwl n 

by (178) and (181). Therefore by (176) and (179), the maximal monotonicity of«H> j3(u) in L p ( R w ) x 
L p (l w ) yields ft = (3(u) in R^ x (0, T). Moreover, /3 (u) coincides with w on l w x (0,T), as w(-,t) 
vanishes outside fl. The rest of proof runs as in the proof of Theorem 2. 


5.5. Proof of Theorem 4. Let s k £ (0,1), Sk \ 0 and let ( u k ,Wk ) be the family of weak solutions 
to the problem 

(182) d t u k + 21 Sk Wk = 0 in X' Skfl , 

(183) w k = %crU k + B(uk) ~ u k in£' a , 


with Ufc(0) = uo,k- Compared to the proofs of Theorems 2 and 3, this proof is definitely easier. 
Actually, since a is kept fixed, the sequence u k retains some space compactness. 

Put X = £ G and suppose that s k < cr for all k £ N without any loss of generality. Thanks to 
the uniform estimates (151)-(153) obtained in §5.1 along with (93), one has, up to a (non-relabeled) 
subsequence, 


Uk —^ U 


21(7 U/s ^ 21(7 u 
Wk —y W 

21 Sk w k -4 W 
/3(u k ) —> P 


weakly star in L°° (0, T; £„), 
weakly in W 1 ’ 2 ^, T; X'), 
weakly star in L°°(0,T;X^ o ), 
weakly in L 2 (0,T; Ho), 
weakly in L 2 (0,T; X'), 
weakly in L 2 (0,T;Ho)> 


which immediately gives w = 21,7 it + /3 — u in £' a . Here we used Lemma 5.1 with /3 > 0 sufficiently 
large so that X' ■—> X' 3 0 (as in §5.3) to identify the limit of fH Sk w k as w. By [36, Theorem 5], one can 
obtain 


Uk^-u strongly in C([0,T];i? 0 ). 

By assumption, uo.fc —>■ uq strongly in Hq. Hence we obtain u(t) -4 uq strongly in Hq as t \ 0. 
Moreover, by applying Minty’s trick to the maximal monotone operator u H > /?(«(•)) in Ho, one 
concludes that j3 = ft(u). For all ip £ ^(1^) C X and <j) £ C'q°(0,T), one can derive 
P T pT 


/ (d t Uk(t),tp)x<l>(t)dt+ (‘<& Sk Wk(t),ip)x dt = 0. 
Jo Jo 


Passing to the limit as k /* oo, we obtain 


f (d t u(t) + w(t), ip) x <£(*) dt = 0. 
Jo 


Since is dense in X, we conclude that 


d t u + w = 0 in X', 0 <t<T. 


Recalling that X = £ a , we conclude that u solves 

dtU + 21<7M + p{u) — u = 0 in £' a , 0 < t < T. 


6. Stationary states 

In this section we analyze the behavior of stationary states of system (l)-(4) in the coercive 
case p > 2. We will put a particular emphasis on the asymptotic behavior of the stationary states 
when er 0. 

The function u is called a stationary state of (l)-(4) when dtU = 0 for a.e. (x, t) € fix (0, +oo). 
Then by (1) we have 

(—A ) s w = 0 in fl, w = 0 in R w \ f 2, 
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which implies w = 0 by Poincare’s inequality (20). Hence, by (2), u solves the problem 

f (~A) a u + fi(u) - u = 0 in ft, 

(184) \ u = 0 in \ ft. 

First of all, let us provide a weak formulation of (184). 

Definition 6.1. A function u € £ a is called a weak solution of (184), if 


(185) 


C(a) 


(u{x) - u{y )) {ip{x) - <p(y)) 


\N+2a 


/K2W \x - y | 

for all <p € £ a . The weak form (185) can be equivalently rewritten as 

21 a u + B{u) — u = 0 in £' a . 


dxdy+ / /3(u)tp dx — / u<p dx = 0 
Jn J n 


We next prove existence of a solution to (184). To this end, we use the direct method of 
calculus of variations. Recall that the energy functional E CT : £ a —> R is defined by 

(186) E a (v):=^p- ff 1 N+ 2 I da;d l/+ / Pi u )& x ~\ f \v\ 2 dx for v G £ a . 

4 JJ R 2 n \x-y\ JS + 2a J a 2Jn 

It is easy to check that E CT is coercive in £ a . Indeed, by Holder’s and Young’s inequalities, 

E„(v) > ^\\v\\ 2 x<r 0 +-^ Jjv\ p dx - C for all v <£ £ a . 

Hence, the existence of a (global) minimizer u £ £ a follows from the compactness of the embedding 
X a Q Hq. Note that u actually solves equation (184) in the sense of Definition 6.1. Now, let u be a 
global minimizer of E^. Since ||u(x)| — |u(y)|| < |u(x) — u{y)\ for (x,y) £ R w x R N , |u| has the same 
energy of u , namely E (T (u) = Eo-dul). Hence |u| also minimizes E CT and solves (184). Consequently, 
by applying maximum principle for the fractional Laplacian (see, e.g., [31], [21]) to the nonnegative 
solution |u|, we infer that |u| > 0 or u = 0 in ft. Therefore every minimizer u of E CT turns out to be 
sign-definite or identically equal to zero over ft. Now, we are going to give conditions implying that 
there exist nontrivial solutions to (184). We can prove the following 


Proposition 6.1. Let Ai(cr) be the first eigenvalue of (48) with r replaced by a. Then, if Xi(a) < 1, 
problem (184) admits a nontrivial weak solution. 


Proof. First of all, we claim that 


inf E ct <0 if Ai(ct) < 1. 

f 

'-'a 


Indeed, let v a be the eigenfunction of 2l CT corresponding to the first eigenvalue Ai(cr), normalized with 
respect to the Ffo- norm - Then it is proved that Ai(er) is simple in [33, Prop. 9 (c)]. Moreover, v a is 
Holder continuous up to the boundary by [32, Theorems 1.1 and 1.3] and [35, Proposition 4]. By a 
simple calculation, we have for e > 0 


M^) = t KII ^ 0 + - / 

z ’ V Jn 

e 2 e p f 

= -77^1 (&) + — \ 

2 P Jn 

7 ; (Al(cr) - 1) + 


\ p dx- 


\ p dx- Y 


' da; 


= £ 


=-P—2 


/ l«.| 

Jn 


1 dx 


< 0 , 


provided that £ is chosen so that 


i(Ar (a)-l) + — [ Krdx<0, i.e. e p ~ 2 < p(1 . 

2 P Jn 2 Jn \v<j\ p dx 

Hence the infimum of E^ over £ a is negative whenever Ai(cr) < 1. Therefore every global minimizer 
is nontrivial. □ 


Remark 6.1. If the first eigenvalue Ai of — A equipped with the homogeneous Dirichlet condition is 
not greater than one, then by the upper estimate Ai(er) < Xf (see (52)), we have Ai(cr) < 1. Hence 
(184) possesses a nontrivial weak solution. 
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Lemma 6.1. Let u be a weak solution of (184). Then its energy is nonpositive, i.e. Eo-(u) < 0. 
Moreover, ifE, a (u) = 0, then u = 0. 

Proof. Let u be a weak solution of (184). Test (184) by u (i.e., substitute ip = u in (185)) to get 

IMlLo+ f \u\ p dx- [ \u\ 2 dx = 0, 

Jn J n 

which yields 

E a (u) = -(^-^j Jju\ p dx<0. 

I 11 particular, if E (T (m) = 0, then u = 0. □ 

As a consequence, we have the following criterion for non-existence of nontrivial solutions: 
Corollary 6.1. If Ai(<r) > 1, then (184) admits only the trivial solution. 

Proof. We observe that, by (49), 

= h\u\\x„ 0 + \ [ \u\ p dx-^ [ \u\ 2 dx 
z P Jn z Jn 

(187) > - (Ai(cr) — 1) / |u| 2 dx+- [ \u\ p dx for all u £ £ a . 

2 Jn P Jn 

Hence, if Ai(tr) > 1, then E ct (m) > 0. Therefore, due to Lemma 6.1, problem (184) has no nontrivial 
weak solution. □ 


We are now in position to state a result on the asymptotic behavior of nontrivial weak solutions as 
a \ 0. 

Proposition 6.2. Suppose that Ai(cr) < 1 for all a £ (0,1). Let, for a £ (0,1), u a he a nontrivial 
weak solution of (184). Then u G converges to zero strongly in Hq as a \ 0. 


Proof. Let a converge to 0 along a sequence crfc and denote by Uk the corresponding nontrivial weak 
solution to (184). Let r > 0 and v £ £ ak be such that |H|# 0 = r. Let us first note that 


\v\ 2 dx < |n|0>- 2 )/ p 



Hence, combining this with (187), we infer 


E<t(i>) > y (Ai(o-fc) - 1) + 


r P 

p | ft |( P " 2 )/ 2 


2 V Al ^ 1 + p|H|(p-2)/ 2 J 


> 0 , 


provided that Ai(o-fe) — 1 + 2r p 2 /(p\Ll\( p 2 P 2 ) > 0, which corresponds to r p 2 > {jp/2)\VL\( p 2 )/ 2 (l— 
Ai(o-fc)). Set 

r k := (|l^l (p - 2)/2 (l - A l( (7 fe ))) 1/(P_2) . 

Then one has 


inf {E ct (i;): v £ £ ak , ||u|| ffo > r k } > 0. 
Therefore, by Lemma 6.1, being u k non trivial, we have 

\\uk\\H 0 < r k - 


Since Ai (a k ) -> 1 as 1; / +00 (cf. Proposition 2.2), we conclude that u k —> 0 strongly in H 0 as 
k /* + 00 . □ 
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7. Appendix 

Here we discuss in some more detail the relations between our problem and other nonlocal 
Allen-Cahn or Cahn-Hilliard models. We just comment on the assumptions on the fractional diffusion 
operators, neglecting the differences occurring in the choice of the potential. 

Let us start describing the relations between our problem and the Cahn-Hilliard equation 
analyzed in [4, 14, 20]. Actually, these problems share a common variational structure as the variable 
w is introduced as the first variation (in L 2 ) of some functional: in our case, of £„■ (cf. (66)), while 
the nonlocal operator J[u] of [4, 14, 20] (recall (7)) corresponds to the gradient of 

(188) Ej[u\ = j [[ j(x-y)(u(x)-u(y)) 2 dxdy, 

4 J Jfixn 

as a direct computation shows. 

Apart from the lower regularity of the kernel K r associated to (—A) 7 ', one major difference 
is that the integration domain in (188) is H x H in place of R 2Ar . In a sense, this is similar to what 
happens in the case of the so-called regional Laplacian (cf., e.g., [23]), defined, for smooth functions, 
as (compare with (10)) 

(189) (-A)? eg u(*) := C(r, N) p.v. y ^+l d V- 

In other words, the above position corresponds to assuming that, as x £ fi, the value of u(y) influences 
that of u{x) only for y £ fi. Hence, in fact, no boundary conditions are taken in that case. In our 
setting, instead, also the outer (Dirichlet) value u(y) = 0, y € R w \ f2, carries some influence on the 
value of (— A.) r u(x). 

This difference is also reflected when one looks at the first variation J of Ej (cf. (7)), where 
the function a(-) depends in fact on the variable x € H. Indeed, taking W = 0 and neglecting the 
principal value for simplicity, also in our case it would be possible to write (at least formally) 

{5E a ){u) = au — K a * u, a= / A cr (a; — y) dy. 

J 

However, a, representing the “total mass” of K ai is now independent of x (note, instead, that the 
convolution term is the same in both cases as we assume u be identically 0 out of f l). 

A further difference between the two models is related to the regularity of the kernels. In 
[4, 14, 20], j satisfies some kind of summability property. On the other hand, in the case of the 
fractional Laplacian, the kernel K r is somehow “less than L 1 ”. This is a trivial remark, of course. 
However, some consequences from the point of view of regularity analysis deserve to be discussed. 
Indeed, in the present case the term 2l r ii is less regular than u and this fact implies that the equation 
enjoys some smoothing effect , as happens in the standard parabolic case (i.e., for the usual Laplacian). 
This fact permits us to “embed” compactness and density tools in the Hilbert formulation by working 
in Hilbert triplets like (X r<0 , H 0: X^ 0 ). 

Instead, in the models analyzed in [4, 14, 20], the term J[u] is strictly more regular than 
the function u on which J acts. This means that the PDE system has limited regularization effects 
(actually, this is especially true for Allen-Cahn based models, see, e.g., [22]; in the Cahn-Hilliard case 
the Laplacian in the equation ut = Am corresponding to our (1) partially compensates this). In other 
words, a singularity in the initial datum tends to propagate with time without smoothing out (at most 
it decreases in amplitude like in a dissipative ODE). 

The regional Laplace operator (189) characterizes also the model studied more recently in [2]. 
To be precise, in [2] the standard Laplacian is taken in the analogue of (1), while the regional operator 
(—A)^ eg appears in the analogue of (2). Moreover, no-flux conditions are assumed for w. This, in 
particular, entails a mass-conservation property; namely, the spatial mean value of u is constant with 
respect to time, as one expects to occur in models describing phase separation, on account of the 
underlying physics. Using that (—A)^ eg can be seen as a fractional power of the Neumann Laplacian 
(restricted to the class of zero-mean functions), the authors of [2] can show that, at least for sufficiently 
smooth solutions, also the component u turns out to satisfy a no-flux condition on dLl (though no 
explicit boundary condition is required in the mathematical formulation of their problem). It is also 
worth mentioning that the analysis given in [2] admits a much wider class of potentials W, including 
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in particular singular functions like the “logarithmic potential” mentioned in the Introduction. In 
principle, this would be possible also for our model and we plan to address this issue in a forthcoming 
work. We also observe that, for the model considered in [2], one expects that the solution, at least 
asymptotically in time and for non-singular functions W, could satisfy strong regularization properties. 
This is due to the fact that the regional Laplacian, as a fractional power of the Neumann Laplacian, 
satisfies the property (—A)£ eg o (—A)£ eg = (—A)^, which may allow use of bootstrap regularity 
methods. On the other hand, an analogue property fails for our operators 2l r due to the occurrence of 
the “solid” Dirichlet condition. For this reason, we expect that the analysis of asymptotic regularity 
properties of weak solutions could be more challenging in our case. 
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